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ABSTRACT 



NOAA-6 satellite AVHRR data and AXBT data were collected 
in the Northeast Pacific Ocean in late 1980 as part of the 
Naval Postgraduate School-sponsored Acoustic Storm Transfer 
and Response Experiment which was in turn part of the U.S.- 
Canadian Storm Transfer and Response Experiment (STREX) . Some 
of the problems in transferring AXBT geographical positions to 
satellite images were solved by designing a computer program 
with accuracies of less than 2 pixels. Thermal comparisons 
were made between AXBT, NOAA-6, and GOSSTCOMP data with the 
result that NOAA-6 data was on the average 2.9°C colder than 
AXBT data and 3.2°C colder than GOSSTCOMP data. Linear regression 
methods reduced to 0.3°C the difference between NOAA-6 and AXBT 
data. Use of this method over a period of 15 days produced a 
mean error of 0.5°C. 

Although NOAA-6 cannot sense directly the subsurface thermal 
structure, it is excellent for observing surface manifestations 
of horizontal thermal features. Further investigation into 
using satellite data as the basis of an empirical relationship 
between the surface temperature and the subsurface vertical 
theinmal structure is warranted. 
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INTRODUCTION 



I . 

The collection of oceanographic data has always been time 
consuming and expensive. With the development of environmental 
satellites, a method for the rapid gathering of oceanographic 
data was available to complement that data gathered on re- 
search cruises. Unfortunately, there are still some problems 
in using satellite data, such as the effects of the atmosphere 
on the radiative transmission path and the effect of the geo- 
metric distortion of thermal features found on polar-orbiting 
satellite images as from NOAA-6. 

This thesis is part of a series of on-going studies at the 
Naval Postgraduate School by the Department of Oceanography 
Environmental Acoustic Research Group. The overall goal of 
the Group is to continue the development of those aspects of 
acoustical oceanography that have a significant effect on naval 
tactical applications. In pursuit of this goal, the Group was 
a participant in the joint U .S . -Canadian Storm Transfer and 
Response Experiment (STREX) held in the fall of 1980 in the 
Northeast Pacific Ocean. The Group is particularly interested 
in investigating whether or not satellites can fulfill the 
role presently played by ships and/or air-dropped expendable 
bathythermographs (AXBT) in gathering sea surface temperature 
data for use in forecasting the ocean's thermal structure, 
using procedures similar to those developed by the U.S. Naval 
Oceanographic Office Antisubmarine Warfare Environmental 
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Prediction (ASWEPS) Program in the early 1960 's. If direct 
correlations could be found between satellite-derived sea sur- 
face temperatures and the vertical thermal structure, then a 
rapid method of surveying the world's oceans could result, 
with numerous naval ramifications. 

As a part of the experiment conducted by the Group under 
the title Acoustic Storm Transfer and Response Experiment 
(ASTREX) , this thesis was directed toward the examination of 
some of the problems in using satellites to observe the hori- 
zontal and vertical thermal distributions in the waters of 
the Northeast Pacific Ocean. Of particular interest was the 
problem of locating open-ocean geographic positions (ship, 
AXBT, etc.) on satellite images with as much accuracy as 
possible. A major portion of this thesis is devoted to this 
subject. If comparisons are to be made between satellite- 
obtained thermal values and ground-truth thermal values, then 
an elimination of location errors between the two media makes 
the results that much more significant. The reasons why NOAA- 
6 satellite imagery is distorted and a method to eliminate any 
location errors successfully are presented below. In addi- 
tion, once location accuracy was assured, various meaningful 
comparisons were made between satellite, bathythermograph, 
and GOSSTCOMP (Global Operational Sea Surface Temperature 
Computation) data. The results of these comparisons are 
also presented. 
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II. REMOTE SENSING IN OCEANOGRAPHIC RESEARCH 



In 1870 and 1879 respectively, the authors Edward Everett 
Hale and Jules Verne wrote about placing an artificial satellite 
in orbit. Hale using a huge waterpowered flywheel and Verne 
a gun of sufficient muzzle velocity (Corliss, 1967). Hale 
envisioned the satellite as an aid to both navigation and 
communication. Up until 1935, the idea of launching man into 
space remained the ideas of small amateur organizations with 
some notable exceptions, such as the efforts of Robert Goddard 
in Auburn, Massachusetts, in the 1920 's. 

With the stirrings of war in Europe in the late 1930 's, 
the now infamous V-2 rockets were developed and launched from 
Peenemuende by a team of German scientists including Wernher 
von Braun. Captured en masse in 1945 by the Allies, this group 
of scientists and their hardware were transferred to the United 
States. The Army Air Force subsequently commissioned a study 
by Rand Corporation, who reported in 1946 for the first time 
that an earth-orbiting artificial satellite could be used 
scientifically in the fields of meteorology, biology, and com- 
munications (Corliss, 1967). This launched a nine-year effort 
by lobbyists both inside and outside the government culminating 
in President Eisenhower's announcement on 29 July 1955 that 
the U.S. would launch an earth-orbiting scientific satellite 
to investigate the environment. The Soviet Union announced 
similar plans on the following day. October 4, 1957 saw the 
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launch of Sputnik-I by Russia followed on 31 January 1958 by 
the launch of the first U.S. satellite, Explorer-I . Since 
that time, hundreds of scientific satellites have been launched 
in order to investigate topics from numerous fields. One of 
these, NOAA-6, was the satellite whose data were used on this 
project. NOAA-6 is an earth-orbiting, environment-sensing 
spacecraft with applications in meteorology and oceanography. 
This section reviews the use of satellites in oceanographic 
studies, provides a description of the operation of NOAA-6, 
and briefly summarizes the normal oceanographic conditions for 
that region of the Northeast Pacific Ocean where the experi- 
ment took place. 

A. OCEANOGRAPHIC CHARACTERISTICS OF THE PROJECT AREA 

The area chosen for this project encompassed that region 
of the Northeast Pacific Ocean between latitudes 40 N and 50 N 
and between longitudes 126 W and 139 W. See Figure 1. The 
oceanographic conditions of this region have been extensively 
studied by Tully (1961; 1964), Tabata (1964; 1965; 1978), and 
Roden (1975) among others. The reader is referred to these 
works for more detailed information as only a brief summary 
of their findings is described below. 

As seen in Figure 1, the project area is located mainly 
in an oceanic water mass transitional region between the Sub- 
arctic Water Mass, predominantly to the north of 45 N, and 
the Pacific Equatorial Water Mass, predominantly to the south 
of 23 N. The Subarctic Current flowing eastward along 45 N 
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Figure 1. Major surface currents and water masses of 

the Northeast Pacific Ocean (after Kibblewhite 
et al. , 1977) 
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latitude divides on the northwest side of the project area 
into the Alaska Current, which circles counterclockwise to 
the north along the Canadian and Alaska coastlines; and into 
the California Current, which flows southward along the western 
coast of the United States. Observations from the project area 
would be expected to show the physical characteristics of the 
Subarctic Water Mass; however, those observations in the 
southern portion of the project area may be somewhat tempered 
by the colder offshore waters of the California Current. 

1 . Thermal-Salinity Structure 

An excess of precipitation over evaporation of approxi- 
mately 25 cm/year (Tabata and Giovando, 1962) has helped to 
create a layer of water extending to a depth of 100 meters in 
the Subarctic Water Mass that is isohaline during the winter 
months. A permanent halocline extending from 100 to between 
200 and 300 meters exists in which the salinity increases by 
1 %o to approximately 33.8 and marks the maximum limit 

of seasonal effects (Tully, 1964) . 

The top of the permanent halocline at 100 meters also 
marks the maximum depth of the seasonal thermocline. During 
the summer, the thermocline forms between 25 and 50 meters, 
influenced heavily by wind mixing effects alone. With the 
coming of the fall and winter months with their intensive 
storms, the surface waters begin to cool considerably and both 
convection and wind mixing erode the thermocline until iso- 
thermal conditions exist to the top of the permanent halocline. 
This condition usually is reached by February at which point 
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the waters continue to cool until the end of March, when the 
heating season begins. See Figure 2 for a general depiction 
of the seasonal structure. Figures 3 and 4 show the expected 
mean thermal structure for the project area for the months of 
November and December. Figure 5 shows the annual surface 
salinity maxima while Figures 6 and 7 show the expected layer 
depths also for November and December. 

Both the Subarctic Current and the California Current 
have surface speeds less than one knot with volume transport 
averaging 10 to 15 million cubic meters per second (Knauss, 

1978) . As a result, the water in the project area is exposed 
to constant climatic conditions over many months and has suffi- 
cient time to adjust to seasonal variations. During November 

and December, there is an expected net heat loss in the pro- 

2 

ject area ranging up to 400 g-cal/cm /day (Tabata, 1961). This 
heat loss is directly responsible for the winter convective 
mixing process mentioned above; therefore, the vertical thermal 
structure is due more to the area's heat budget, storm cycle, 
and salinity layers than to any influx of new water. 

2 . Internal Waves 

Internal waves have been shown to cause amplitude 
fluctuations of 4.5 to 5 meters at the level of the seasonal 
thermocline in the Northeast Pacific Ocean (Tabata and Giovando, 
1962) . These oscillations may cause a periodic thickening of 
the thermocline from 10 to 50 meters due to phase differentials 
between the top and bottom of the thermocline (Tully, 1964) . 
These periodic fluctuations vary with depth as oscillations of 
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Figure 2 . Seasonal oceanographic structure of the 
Northeast Pacific Ocean (from Tully, 
1964 ) 
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Figure 5. Annual surface salinity of the project 
area (from Robinson, 1976) 
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Figure 6. November mean layer depth in project 
area (from Robinson, 1976) 
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Figure 7. December mean layer depth in project 
area (from Robinson, 1976) 
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isotherms at the level of the halocline are five times 
greater than oscillations in the region of the thermocline. 

3 . Fronts and Eddies 

The Subarctic Front, usually found between latitudes 
40 N and 45 N, may be present in the center of this project 
area. This front is characterized by the lack of a density 
front in the upper 100 meters, by the region of the strongest 
surface baroclinic flow being to the south of the surface 
temperature and salinity fronts, and by the mixed layer depth 
extending to the top of the halocline at 100 meters on its 
northern edge (Roden, 197 5) . In the southern section of the 
project area, eddy formation may be present in that area in- 
fluenced by the California Current which has been described 
as a relatively shallow meandering current with alternating 
warm and cold tongues (Bernstein ^ , 1977). 

B. USE OF SATELLITES IN OCEAN THERMAL STUDIES 

In 1968, a study was done comparing satellite-obtained sea 
surface temperatures with monthly mean surface temperatures 
with the result that the satellite values were anywhere from 
3 to 8.3 degrees C lower than the mean (LaViolette and Chabot, 
1968) . The relative horizontal gradients observed in their 
satellite data, however, were fairly consistent with similar 
mean gradients from the historical data. This pattern of 
satellite-obtained sea surface temperatures being lower than 
the mean or the actually observed sea surface temperatures 
persists until today, except that technological advances have 
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reduced the differences in temperature between the two sets 
of data so they now are between 0.5 to 3.0 degrees C (Rao et 
al . , 1972; Brower et , 1975; McMillin, 1975; Cogan and 
Willard, 1976; Barnett e^ , 1977; Tabata and Gower, 1980). 

Prior to 1972, the oceanographic use of satellite-obtained 
sea surface temperature was severely limited by both the en- 
gineering characteristics of the satellite radiometers and by 
the environmental aspects causing atmospheric attenuation. 
Large instantaneous f ields-of-view (IFOV) limited the resolu- 
tion capacity of the satellite and the large values for the 
variations in the electronic signal (NEAT) caused spatial and 
temporal errors, making it difficult to detect the gradients 
associated with oceanic fronts (Legeckis, 1978) . 

Several methods were proposed to remove the atmospheric 
contamination responsible for the majority of the difference 
between satellite and observed sea surface temperature values. 
The 3 to 8.3 degree difference found by LaViolette and Chabot 
(1968) came from satellite data that were not corrected for 
atmospheric attenuation, but in 1969 they developed a daily 
averaging method to lessen its impact (LaViolette and Chabot, 
1969). Vukovich (1971) developed a filtering technique to 
accomplish the same purpose while Smith e^ al . , (1970) used 

a statistical method which, when compared with ship observa- 
tions, had both bias and random errors of less than 1 degree 
C using early NIMBUS satellite data. Maul and Sidran (1973) 
investigated the effects of the atmosphere, nadir angle, cloud 
amount, cloud height, and random noise which resulted in a 
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theoretical error (2 degrees C) for the NOAA satellite series, 
then soon to be launched. 

In early 1970, NOAA launched ITOS-1 which was the first 
satellite in the NOAA series of satellites of which both NOAA- 
6 and NOAA-7 are now in orbit. In early 1970, NOAA-NESS be- 
gan working on a satellite data processing model, to include 
the effects of atmospheric attenuation, which was the prede- 
cessor to the GOSSTCOMP (Global Operational Sea Surface Tempera- 
ture Computation) model (Brower et a^. , 197 6). With the 
launch of NOAA-2 in 1972, a more advanced radiometer was put 
into use with an IFOV of about 1 kilometer and a much reduced 
system NEAT of less than 3.0 degrees C (Legeckis, 1978). With 
this improved system, sea surface temperature fronts could be 
detected and monitored. Among the studies done during the 
following few years were those of LaViolette (1974) on upwell- 
ings off the west coast of Africa, Stumpf and Rao (1975) on 
tracking eddies in the Gulf Stream, and Bernstein e^ a_l . , (1977) 

on the comparison of eddies in the California Current with 
direct observations. 

By the mid-1970's NOAA-NESS had refined their satellite 
data processing model; however, comparisons with observed data 
by NOAA itself and by others found that the quality of measure- 
ments varied with time and geographical area and were related 
to the temperature gradient field; good correlation came from 
regions of weak gradients and marginal results came from regions 
of strong gradients (Brower e^ a]^. , 1976) . Klein (1979) found 
that NOAA-5 sea surface temperatures in the Northeast Pacific 
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Ocean that had been subjected to the GOSSTCOMP model were 
biased 3.5 to 3.9 degrees C and suggested that the error was 
a result of overcorrection by the model for atmospheric attenu- 
ation. With the launching of TIROS-N in 1978, NOAA-NESS up- 
dated GOSSTCOMP to take advantage of the Advanced Very High 
Resolution Radiometer (AVHRR) on this, and on the follow-on 
NOAA-6 and NOAA-7, satellites. Whereas NOAA-5 had a NEAT of 
1 to 1.5 degrees C, the TIROS-N/NOAA A-G satellite series has 
a NEAT of 0.12 degrees C (Schwalb, 1978). The improvement in 
NEAT should result in a better correlation between observed 
and satellite-derived sea surface temperatures. Chahine (1980) 
suggested that an absolute accuracy of 1 degree C in these 
differences could be obtained by simultaneous observations of 
atmospheric and surface emissions with multi-channel radiometers, 
using spectral regions of the 3.7 vm carbon dioxide windows as 
the main sounding channel. An instrument to accomplish this 
has yet to fly on a satellite. 

1 . Problems Associated with a Satellite Data Base 

Briefly described below are three common problems asso- 
ciated with using satellites in thermal studies. Atmospheric 
attenuation is important when interpreting satellite-derived 
temperatures, while location accuracy is important when ther- 
mal comparisons are made between ship, satellite, and AXBT 
data. The depth to which present-day radiometers sense the 
thermal structure concludes the section. 
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a. Atmospheric Attenuation 

As will be described in detail in Section II. C. 2, 
data from the 10.5 to 11.5 pm infrared channel on NOAA-6 were 
used on this project. Radiation in this spectral region 
emitted from the earth's surface or from cloud tops is attenu- 
ated in its passage through the atmosphere to the radiometer. 

The major contribution to this attenuation is water vapor which 
can be responsible for up to a 9.0 degree C correction in the 
satellite data (Brower al^. , 1976) . The amo\int of water 
vapor in the atmosphere varies horizontally, vertically, and 
in time with the least amount of absorption around the 9.5 to 
10.5 pm region (Fett and Mitchell, 1977). Other absorbers and 
their possible corrections are carbon dioxide (0.1 to 0.2 
degrees), ozone (0.1 degrees) , and aerosols (0.1 to 0.95 degrees) . 
Details on the physics of this absorption process can be found 
in Roberts ^ ^. , (1976) and Weinreb and Hill (1980). 

Many atmospheric correction techniques have been 
tried in an attempt to correct satellite data. Some of these 
were discussed previously. A knowledge of the vertical mois- 
ture field would help significantly in reducing the attenuation 
effects but these data are not generally available. In any 
case, the multispectral approach to this problem seems to offer 
the best chance to reduce this type of error significantly 
(Chahine, 1980; Deschamps and Phulpin, 1980). 

b. Location Accuracy 

A major portion of this project was devoted to 
locating geographic positions correctly on satellite imagery. 
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Estimates of location error vary widely. The data archived 
from early satellites in the NOAA series had accuracies within 
20 kilometers along the orbital track but for high zenith 
angles the accuracy decreased to within 40 kilometers (Conlon, 
1973) . Subpixel accuracies were discussed by Bernstein and 
Ferneyhough (1975) on LANDSAT imagery. A satellite image con- 
taining land is usually much easier to correct geometrically 
than an image that contains open ocean, simply because control 
points are much easier to identify on the land image. 

A method of correcting VHRR imagery using a simple 
algorithm by Legeckis and Pritchard (1976) had a mean accuracy 
of 5 kilometers. A technique whereby ship positions were 
transferred to satellite images by comparing cloud features 
yielded errors of 80 to 90 kilometers although this was not the 
main purpose of the study (Cogan and Willand, 1976) . Another 
study, using both a zoom transfer scope and NOAA-supplied photo- 
graphic prints of satellite imagery used in conjunction with 
plastic overlays containing latitude and longitude lines, had 
a location accuracy to within 10 to 50 kilometers depending on 
the distance of the feature being located from land control 
points (Tabata and Gower, 1980) . 

A full discussion of location accuracy can be found 
in Section IV below. 

c. Remote Sensing of the Vertical Structure 

The sea "surface" temperature that a satellite 
senses is a manifestation of a complex process that occurs in 
the top few millimeters of water called the thermal boundary 
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layer. This layer is subject to the processes of net upward 
heat flux, infrared and solar radiation, and turbulence with 
the resulting temperature difference between the top and bottom 
of this layer of up to 1.0 degrees C (Katsaros, 1980). Typical 
radiometers sense only the radiation emitted from a depth of 
about 50 vm. 

Direct measurement by satellites of the deeper 
vertical thermal structure is not possible with the instriiments 
carried onboard the satellites in orbit today. Techniques 
using Raman lidar systems have been developed theoretically 
and prototypes experimentally tested with reported accuracies 
within 0.2 degrees C (theoretical best value) to depths of 30 
meters (Leonard e^ a_l. , 1979) . Conclusions from this study 
suggest that the structure to depths of 100 meters may be de- 
tectable. The physics of the Raman spectra used in this pro- 
cess can be found in Murphy and Bernstein (1972) . 

C. NOAA-6 OPERATION 

The NOAA-6 satellite is the second satellite in a series 
of third generation, polar-orbiting satellites that began with 
the launch of TIROS-N on 13 October 1978 at the Air Force 
Western Test Range, Vandenberg Air Force Base, California. The 
TIROS-N/NOAA A-G satellite series, of which NOAA-A became re- 
designated NOAA-6 upon its successful launch, is a joint re- 
search effort of the United States, the United Kingdom, and 
France and is operated by the National Environmental Satellite 
Service of the National Oceanic and Atmospheric Administration 
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(NOAA-NESS) under the U.S. Department of Commerce. The United 
Kingdom provided one of the three sounding units onboard the 
satellite, France supplied the onboard data collection system 
(DCS), the National Aeronautics and Space Administration (NASA) 
funded the development and launch of TIROS-N, and NOAA supplied 
the funds for the NOAA-6 satellite. The mission objective of 
this satellite series that directly relates to this thesis is 
the continuous monitoring of the environmental features in the 
western hemisphere which is accomplished in conjunction with 
a second satellite system, also operated by NOAA, the Geo- 
stationary Operational Environmental Satellite (GOES) System. 

It should be noted that TIROS-N ceased operation in late 1980. 
NOAA-6 was still functioning at the writing of this thesis and 
NOAA-7 began operating in June 1981. 

For the purposes of this project, only those spacecraft 
systems that were extensively used or are important to the 
understanding of the results are explained below. The reader 
is referred to Schwalb (1978), Hussey (1979), Lauritson, et al . , 
(1979) , and ITT Aerospace (undated) for a fully detailed des- 
cription of the many instruments onboard NOAA-6. Sections of 
these references, especially the works of Schwalb and Hussey, 
were used extensively below. 

1. The Spacecraft 

NOAA-6 used an Atlas-F launch vehicle which is a com- 
paratively small rocket approximately 28 meters tall and 
weighing about 600,000 kilograms. See Figure 8. The main body 
of the rocket detaches after launch and a second stage solid 
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Figure 8. Atlas-F launch vehicle (from Hussey, 1979) 
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rocket motor, an integral part of the NOAA-6 satellite itself, 
burns until depletion putting the satellite into a nominal 
833-kilometer orbit. 

a. Physical Structure 

The satellite itself, as shown in Figure 9, con- 
sists of three sections. The Reaction Support Structure (RSS) 
includes the injection motor mentioned above, the attitude 
control propulsion system, and an 11.6 square-meter solar cell 
array. The Instrument Mounting Platform (IMP) includes the 
attitude control sensors and the Advanced Very High Resolution 
Radiometer (AVHRR) . The five-sided central structure, located 
between the RSS and the IMP, includes twelve thermal control 
louvres and the earth-facing communications antennae. The 
satellite is 3.71 meters long and 1.88 meters in diameter. Its 
weight at laxinch was 1420 kilograms which reduced to 7 37 kilo- 
grams once established in its orbit. 

b. The Attitude Determination and Control Subsystem 

(ABACS) 

When a satellite sensor, such as the AVHRR, scans 
the surface of the earth, the attitude of the spacecraft is 
extremely important in determining during data analysis just 
where the sensor looked. Any roll, pitch, or yaw on the sat- 
ellite will make the application of scan geometry extremely 
difficult and significant errors would result. Because of 
this, the ABACS system was designed to maintain the attitude 
of the spacecraft to within 0.2 degrees (3-sigma) of the local 
geographic reference (Schwalb, 1978). This value is obtained 
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through the use of three mutually-orthogonal torque wheels 
which receive input from the Earth Sensor Assembly (ESA) for 
pitch and roll and, for yaw, an inertial reference source with 
sun-sensor updates. The ESA is an infrared (IR) sensor that 
views the entire earth and supplies torque input to keep the 
earth centered between four independent detectors. The sun 
sensor uses multiple data inputs from various mechanisms to 
provide the yaw input. 

c. Data Handling Subsystem 

There are four primary components in the data 
handling system onboard NOAA-6; the TIROS Information Processor 
(TIP) , the Manipulated Information Rate Processor (MIRP) , the 
Digital Tape Recorders (DTR) , and the Cross Strap Unit (XSU) . 
All the information eventually received on the ground from 
NOAA-6 has to be processed by at least one of these four com- 
ponents. Figure 10 is the data flow diagram for NOAA-6; atten- 
tion is drawn to the path followed by the AVHRR data via the 
MIRP to the switching unit for subsequent transmission at 0.66 
megabits per second to the earth station antenna as real-time 
High Resolution Picture Transmission (HRPT) data. The AVHRR- 
HRPT data stream was the only one used on this project. An 
explanation of the various other sensors shown in the diagram 
can be found in Schwalb (1978) . 

The MIRP formats the AVHRR data and adds synchroni- 
zation, identification, telemetry, and time code information. 

It senses a pulse at the start of each AVHRR scan line and 
initiates a data sampling process that divides the arriving 
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1702.5 MHf. LEFT HAND CIRCULAR 
AVAILABLE BUT NOT PLANNED 
FOR HRPT USE 



earth scan data into 2048 computer data words per scan line. 

The pulse that is sensed at the initiation of each scan line 
originates when the AVHRR scan mirror, which rotates at 360 
RPM producing 6 scan lines per second, reaches a precise posi- 
tion in its sweep just prior to scanning across the surface of 
the earth. The data are stored in memory and then subsequently 
read out at a rate suitable for the HRPT on a first-in first- 
out basis. Any one of the 2048 data words or samples, along 
with the number of the scan line on which it is located, de- 
fines a pixel. Throughout this project, the term pixel will 
be defined by the designation (scan line number , sample number) 
or, in short, (NL,NS) . A more comprehensive discussion of this 
process can be found in Section III.B below. 

2 . NOAA-6 Onboard Sensors 

There are three primary environmental sensors onboard 
NOAA-6. The TIROS Operational Vertical Sounder (TOVS) con- 
sists of the High Resolution Infrared Radiation Sounder (HIRS/ 
2) whose purpose is to provide data to allow calculation of 
vertical temperature profiles and atmospheric water and ozone 
concentrations, the Stratospheric Sounding Unit (SSU) , and the 
Microwave Sounding Unit (MSU) . The second sensor, the Space 
Environment Monitor (SEM) , consists of a Total Energy Detector 
(TED) , the Medium Energy Proton and Electron Detector (MEPED) , 
and the High Energy Proton and Alpha Detector (HEPAD) . The 
last of the three systems, the AVHRR, was the sensor system 
extensively used on this project and will be described in 
detail below. For an in-depth discussion of the first two 
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sensor systems mentioned above, the reader again is referred 
to Schwalb (1978) . 

The AVHRR aboard NOAA-6 is a four-channel scanning 
radiometer that is sensitive to energy in four regions of 
the electromagnetic spectrum. Table 1 below is a summary 
of NOAA-6 channelization. 



Table 1 

NOAA-6 AVHRR Channelization 
(adapted from Schwalb, 1978) 



CHANNEL 


WAVELENGTH (um) 


REGION 


PURPOSE 


1 


0.58 - 


0.68 


visible 


cloud coverage 
land-water bound, 
snow-ice extent 


2 


0.725 


- 1.1 


visible- 

near-ir 


as above 


3* 


3.55 - 


3.93 


mid-ir 


sea surface temp, 
cloud mapping 


4 

* 


10.5 - 


11.5 


far-ir 


sea surface temp, 
cloud mapping 



On NOAA-6, channel 3 is very noisy and usually not used 



An afocal 20.3 cm-aperture telescope, which produces 
a field of view of 1.3 ± 0.1 milliradians and an instantaneous 
field of view (IFOV) ground resolution of 1.1 kilometers at 
nadir (Lauritson, ^ a^. , 1979) , separates the radiant energy 
into the four spectral regions with the help of secondary op- 
tics. The radiant energy in each of these regions is then 
focused on its respective detector. The quantity of energy 
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sensed by the detector then is converted to a count value from 0 
to 255 in the format used for this project. Channel 4 was the 
main channel from which information was gathered and it uses 
a mercury cadmium teluride (H^C^T^) detector optimized for 
best sensitivity between 10.5 and 11.4 micrometers (Schwalb, 

1978) . The spectral response curve for channel 4 is shown in 
Figure 11. The noise equivalent differential temperature 
(NEAT) , a measure of the random or coherent two-dimensional 
noise patterns superimposed on the data signal broadcast to 
earth, is less than 0.12 degrees Kelvin at 300 degrees Kelvin. 

Pre-launch AVHRR calibration is covered in a report 
by ITT Aerospace (undated) and post-launch thermal calibration 
of channel 4 is covered extensively in a report by Laurtison, 
et al. , (1979). For every scan line, the radiometer views deep 

space (0 radiance) and then a blackbody target designed into 
the radiometer housing and kept heated to 15 degrees centigrade. 
To a first order approximation, the radiometer output is linear 
with input energy (Schwalb, 1978) so a two-point linear cali- 
bration, using the above values, is done during every scan 
sequence. Channel 4 with its H^C^T^ detector, however, has a 
not-quite-linear response due to the physical properties of 
the HgC^T^. Lauritson, et ^. , (1979) have generated a table 

of errors for this channel, which represents the difference 
between the actual target temperature and the temperature de- 
rived from the two-point calibration. Table 2 is a summary of 
these data. Note particularly the small errors around 285 
degrees Kelvin, for this is the sea surface temperature range 
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11. NOAA-6 AVHRR channel 4 spectral response 
curve (from Kidwell, 1979) 
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determined by the AXBT drops. With this information, a table 
of count-value-to-temperature conversions was generated for 
the time of this project by NOAA-NESS and is included in 
Appendix A. 



Table 2 

NOAA-6 AVHRR Channel 4 Nonlinearity Errors 
(from Lauritson e^ , 1979) 

TARGET TEMPERATURE ERROR 

(degrees K) 



305 


0.5 


295 


0.3 


285 


O 

• 

o 


275 


• 

o 

1 


265 


1 

o 

• 

CO 



3. NOAA-6 Orbital Parameters 

Because the development of the satellite data set was 
so closely intertwined with the NOAA-6 orbital parameters, 
their discussion is included in Section III.B.3.C below. 
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Ill . DATA COLLECTION AND PROCESSING TECHNIQUES 



The use of satellite data on any research topic intro- 
duces extensive data processing problems, especially when one 
considers that a typical NOAA-6 infrared satellite image con- 
tains over nine million pieces of data. This section explains 
the procedures used on this project to collect, process, and 
analyze NOAA-6 satellite imagery with emphasis in the area 
of geographic location accuracy. Also included in this sec- 
tion are the procedures to collect and process the AXBT data 
as well as the collection of the GOSSTCOMP product. 

A. AXBT COLLECTION AND PROCESSING 

As part of the data base for this project, a series of six 
Navy P-3C aircraft flights were staged out of NAS Moffett Field, 
California, for the purpose of dropping a pattern of bathyther- 
mographic sonobuoys (AXBT) . The dates of these flights were 
15, 17, and 19 November and 1, 3, and 5 December 1980. These 
flights were scheduled as part of the Naval Postgraduate School's 
research effort on behalf of the joint U . S . -Canadian Storm 
Transfer and Response Experiment (STREX) . 

Eact of the nine-hour flights flew northwestward from Cape 
Mendocino, California and proceeded to drop a series of AXBT's 
along a track 1333 kilometers (720 nm) long as shown in Figure 
12. The spacing between the buoys was 55.6 kilometers (30 nm) . 
The first four flights flew out and back on the center track 
dropping AXBT's on positions 1 through 24. Flight 5 flew the 
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Figure 12. AXBT patterns for the project area 
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center track on the outbound leg and flew the southern track 
on the return leg dropping AXBT's on positions 1 through 13 
and 34 through 46. Flight 6 repeated the center track outbound 
and flew the northern track inbound dropping AXBT's on positions 
1 through 13 and 54 through 65. The northern and southern 
return tracks were designed to gather data on the horizontal 
thermal structure and were offset 111 kilometers (60 nm) 
either side paralleling the center track. 

The complete navigation suite of the P-3C was used in cal- 
culating the position of each of the deployed AXBT's. At the 
end of each flight, the cumulative error of the inertial navi- 
gation system was checked and recorded. For the flights whose 
data were selected for the project, this error was less than 
4 nautical miles . 

Also important to note is the procedure that the P-3C on- 
board computer uses to calculate the splash points of the de- 
ployed AXBT's. Upon releasing the AXBT from the aircraft, the 
computer ballistics program uses the calculated wind speed 
and wind direction from the navigation system to provide a 
trajectory for the first 2000 feet of fall. After this 2000 
feet of fall, the ballistics program assumes a straight descent 
to the water. This entry point becomes the so-called splash 
point for which geographical coordinates are calculated and 
displayed to the flight crews. Most of the AXBT's were dropped 
from an altitude of 2000 feet except when low clouds or icing 
conditions prevented flying at that altitude. It was felt 
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that the location error associated with the few high altitude 
drops was within the 4 nm aircraft navigation error. 

1 . The Air-Dropped Expendable Bathythermograph (AXBT) 

The AXBT is an air-dropped expendable bathythermo- 
graph transmitter set deployed by Navy P-3C and S-3A aircraft. 
Its purpose is to provide an accurate profile of the vertical 
thermal structure from the ocean's "surface" to about 350 
meters. Upon water entry, a seawater battery activates, power- 
ing a VHF transmitter, and approximately 30 seconds later a 
temperature probe begins a 5 foot/second descent (Sparton 
Electronics, 1976). The temperature probe and accompanying 
electronics within the sonobuoy package translate the sensed 
water temperature into a frequency broadcast by the radio trans- 
mitter using the formula 

frequency = 800 + 20 (temperature deg. F.). 

This low-power broadcast from the sonobuoy is intercepted by 
the aircraft where it is electronically recorded on specially 
processed paper in real-time. 

The accuracy of this process is governed by the accuracy 
of the thermal probe on the sonobuoy and this is claimed to be 
within 1.0 degrees C by the manufacturer (Sparton Electronics, 
1976) . Reports in the literature place the repeatable accuracy 
to within 0.2 degrees (Barnett ^. , 1979). The detailed 
workings of this type of AXBT is described in Sessions and 
Wilson (1976) although the buoy described in their work was 
supplied by a different manufacturer. 
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It is also important to note that the temperature probe 
does not start at the surface of the ocean but begins its 
descent from a depth of about 0.2 meters (Barnett ejt a]^. , 1979). 
The probe itself also may be subject to very low temperatures 
if the aircraft transports the buoys at a high altitude for 
a long period of time before deployment or if the buoys them- 
selves are dropped from a high altitude. Both Navy aircraft 
described above have systems designed to prevent the freeze-up 
of the sonobuoys, and altitude launch restrictions do apply 
for the deployment of this buoy. 

2 . AXBT Data Processing 

The thermal profiles were recorded using two different 
methods. As described above, the P-3C-produced paper copy of 
the thermal profile was used with a plastic overlay to read 
off the temperature for any depth. These readings were then 
transferred to paper logs by hand. The accuracy of this method 
is within the accuracy limits of the AXBT itself. 

The second method involved the use of an AXBT-digitizer 
provided by the University of Hawaii and which was used also 
in the NORPAX Experiments. This piece of equipment was de- 
signed specifically to be used onboard the P-3C aircraft during 
flights. It connects into the equipment that receives the 
signal broadcast from the AXBT and digitally records on mag- 
netic tape the signal representing the thermal profile at one 
second intervals. These tapes then are analyzed on a computer 
and various outputs produced. Figure 13 is an example of a 
group of AXBT profiles produced by this system. Figures 14 and 
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Figure 13. Digitizer AXBT profiles, an example 
(from Kilonsky, 1981) 
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Figure 14. Digitizer depth-of-isotherra summary, an example 
(from Kilonsky, 1981) 
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15 show two additional products, a depth-of-isotherm summary 
and a listing of the temperature at 5 meter intervals for 
each AXBT. The accuracy of this system is also within the 
accuracy limits of the AXBT. 

It should be noted also that the AXBT-digitizer had 
provisions for recording the raw received signal from the 
AXBT onto an analog tape recorder. These tapes then were 
used during the analysis phase as a direct input to the 
AXBT-digitizer in order to verify questionable temperature 
profiles . 

B. SATELLITE DATA SET SELECTION AND PROCESSING 

The decision criteria used to determine which satellite 
passes to examine were reviewed in the following order: 

(1) the satellite pass coverage had to include the ocean 
area where the AXBT's were dropped; 

(2) the time of the satellite pass should be as close as 
possible to the time when the AXBT's were dropped; 

(3) the ocean areas containing the AXBT's should be rela- 
tively cloud-free; and 

(4) there had to be at least one clearly identifiable 
landmark somewhere on the full satellite image. 

Two outside government facilities were used in addition to 
the facilities at the Naval Postgraduate School in order to 
choose satellite passes which met these decision criteria. 

1. NOAA-NESS 

The facilities of the Satellite Field Services Station 
of the National Oceanic and Atmospheric Administration's National 
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Figure 15 



Digitizer temperature- at -5 -meter-intervals 
summary, an example (from Kilon'sky, 1981)' 



date 




time 


iat. long. 




axbt 


depth 


temp 


1. 






\ 




\ 










no. 

1 


(m) 


(C) 








1191180 


3 


192014 


425500 1291200 


50 


1 

4 


1 


1 






77 


01481 


514IT3 101483 


151403 


201401 


231403 


"301401 


33T402 


401402 


451482: 


50K02 


551421 601341 


651217 


701 104 


731071 


801045 


831032 


901013 


95 972 


100 9o3 


105 


945 110 


929 


115 917 


120 


91 I 


123 904 


130 902 


135 889 


140 


881 


145 062 


150 831 


1 55 


813 lOO 


030 


105 033 


170 


827 


175 030 


180 819 


185 81 1 


190 


797 


195 794 


200 703 


203 


780 210 


770 


215 773 


220 


770 


225 705 


239 700 


233 752 


240 


740 


245 740 


250 733 


255 


730 200 


725 


205 717 


270 


71 I 


273 705 


289 098 


283 094 


290 


088 


295 682 


300 ’<>74 


305 


07 1 3 1 0 


003 


315 002 


320 


037 


325 055 


331) 052 


335 049 


340 


644 


345 644 


350 030 


335 


030 iOO 


02 V 


365 022 


370 




375 014 


389 612 










2191130 


3 


I 933111 


433800 130 


800 


59 


5 










77 


01 V-/3 


51 4'* 5 101443 


151440 


201440 


251430 


30 1 440 


331442 


401440 


451442 


30M46 


351419 001311 


051130 


701091 


731059 


891049 


051034 


901021 


95 996 


100 9C0 


105 


959 110 


9 t3 


1 1 5 920 


120 


912 


125 897 


130 900 


133 894 


140 


880 


145 878 


1 oO oo j 


155 


8'3 100 


329 


105 822 


170 


810 


175 81 I 


181) 803 


185 797 


190 


789 


195 703 


200 772 


205 


701) 210 


770 


215 770 


220 


7 08 


225 702 


230 757 


235 75*1 


240 


751 


243 745 


250 723 


255 


735 200 


727 


205 722 


270 


719 


275 713 


289 708 


280 700 


290 


092 


295 607 


CCO 070 


305 


071 310 


003 


315 052 


320 


044 


325 034 


330 025 


335 014 


340 


009 


345 604 


350 390 


3o 5 


53') 3 60 


507 


305 504 


370 


577 


375 574 380 571 










31^110^ 


3 


1 95034 


412100 131 


700 


55 


0 










77 


01303 


51390 101397 


151395 


201395 


251395 


301397 


351393 


401304 


431394 


501305 


55 1 


222 001120 


031039 


70 


991 


75 907 


80 959 


85 949 


90 


932 


93 918 


100 003 


105 


9C0 1 1 0 


370 


1 1 5 8T0 


120 


802 


123 853 


130 838 


135 827 


140 


814 


145 701 


1 oO 7 3 


155 


705 100 


701 


105 752 


170 


744 


175 738 


180 735 


183 738 


190 


733 


195 734 


200 733 


205 


735 210 


727 


215 727 


220 


727 


225 719 


230 717 


235 702 


240 


094 


245 692 


230 024 


255 


602 200 


071 


205 003 


270 


003 


275 05 I 


280 040 


285 044 


290 


039 


295 636 


200 023 


305 


617 310 


006 


315 000 


320 


590 


325 393 


330 581> 


333 582 


340 


574 


345 571 


330 371 


355 


505 700 


501 


305 555 


370 


550 375 547 


380 545 










IT17 




2J"258 


45> 


400 132 


700 


58 


7 










77 


01314 


5131 1 101314 


151314 


201317 


251319 


301317 


351310 


401316 


451314 


501317 


551 SCO 001 112 


051013 


70 


988 


75 952 


80 933 


85 920 


90 


900 


95 009 


IGG :70 


105 


C5l 1 1 0 


835 


115 814 


120 


794 


125 783 


130 709 


133 754 


149 


744 


145 733 


130 723 


1 55 


733 1 00 


733 


105 733 


170 


733 


l75 735 


180 733 


1 85 < 3 v> 


190 


735 


195 733 


200 724 


205 


732 210 


727 


215 727 


2‘*0 


722 


225 722 


230 722 


233 718 


240 


712 


245 710 


250 701 


255 


095 200 


090 


205 084 


270 


082 


275 077 


280 007 


285 037 


290 


0>7 


295 639 


300 031 


305 


023 310 


020 


315 013 


320 


008 


325 001 


330 590 


335 582 


340 


574 


343 371 


330 703 


355 


557 300 


550 


303 547 


370 


:;>!3 


375 536 


380 536 










ol6l IGO 




20 1 523 


4540CO 133 


800 


03 


0 










77 


012C0 


51300 IOI300 


151390 


201300 


251303 


301303 


351300 


401305 


451306 


3012G3 


551300 001 


237 


051029 


70 


904 


75 922 


80 900 


83 885 


90 


870 


95 039 


100 234 


105 


IH J 1 1 0 


827 


I 1 3 790 


120 


709 


125 742 


130 719 


135 793 


140 


098 


145 098 


130 093 


1 c ^ 


701 100 


702 


105 702 


170 


098 


175 702 


180 705 


185 700 


100 


706 


195 703 


200 702 


205 


701 210 


095 


215 095 


220 


090 


225 087 


230 084 


235 679 


240 


074 


245 668 


250 003 


255 


000 200 


052 


205 019 


270 


049 


275 041 


280 035 


283 031 


290 


928 


293 625 


3C0 020 


3C5 


012 310 


012 


315 001 


320 


395 


325 599 


330 5U2 


335 574 


340 


571 


345 369 


330 TOO 


354 


557 ‘HO. 


553 


305__5:?0. 


370 


51 


375 539_ 


380__531 










6 1 9 1 I 20 


3 


292220 


402300 1341200 


00 


9 










77 


01 193 


51 


203 1 0 I 


204 


1 5 1 204 


201 


20 !• 


251204 


301204 


351203 


401 


203 


451204 


50 1 204 


551 


20 ’ oO 1 


201 


051 130 


70 


833 


75 050 


80 837 


85 824 


90 


81 1 


93 805 


100 30 


105 


797 110 


734 


I I 5 709 


12') 


75JI 


125 742 


1'30 095 


135 081 


140 


079 


145 684 


130 090 


155 


094 1 oO 


701 


105 701 


170 


703 


175 703 


180 701 


183 703 


190 


701 


195 698 


200 s>95 


205 


08o 210 


079 


215 082 


220 


082 


225 074 


239 071 


235 097 


240 


065 


245 663 


250 OoO 


255 


055 200 


049 


205 041 


270 


(i39 


275 028 


289 023 


285 017 


290 


012 


295 606* 


3G0 001 


305 


590 310 


590 


315 581 


320 


570 


325 509 


330 500 


335 555 


340 


v) 53 


343 347 


310_ 74 5 


355 


53'^ :]<-0 




305 527 


370 




'375 518 


380 512 










7191 ICO 


3 


204123 


46 


700 1334000 


04 


17 










77 


01279 


51202 101 


282 


151282 


201 




25 1 282 


301202 


331282 


401 


201 


451201 


50 1 2'*2 


551 


270 001 


271 


(>5I 147 


701031 


75 948 


80 995 


83 885 


90 


870 


95 857 


100 B-:o 


105 


832 1 10 


819 


1 1 5 000 


120 


788 


123 772 


139 757 


135 751 


140 


744 


145 730 


150 722 


15 5 


721 100 


I 


105 724 


170 


730 


175 727 


180 727 


185 735 


190 


7 '30 


195 719 


200 719 


205 


718 210 


7T9 


213 717 


220 


717 


225 714 


2*30 71 1 


235 701 


240 


097 


245 690 


230 


255 


032 2oO 


071 


205 oOO 


270 


057 


275 057 


280 052 


285 047 


290 


036 


295 031 


300 023 


305 


020 310 


017 


313 009 


320 


005 


325 001 


339 593 


333 588 


340 


583 


345 577 


350 574 


355 


503 300 


504 


305 537 


370 


553 


375 545 


380 542 











56 



Environmental Satellite Service (NOAA-NESS) in Redwood City, 
California, were used in initially selecting the satellite 
passes. The Redwood City facility is one of three NOAA-NESS 
stations that monitor NOAA-6; the other two are the Command 
and Data Acquisition (CDA) stations in Gilmore Creek, Alaska, 
and Wallops Island, Virginia. Redwood City differs from the 
CDA stations in that Redwood City records the digital High 
Resolution Picture Transmission (HRPT) readout consisting of 
three channels of AVHRR data in the 8-bit precision field- 
station format. These 1600 BPI, 9-track computer-compatible 
magnetic tapes then are archived in Redwood City on a 90-day 
rotating basis. The CDA stations record the various other 
data formats broadcast from NOAA-6 as well as recording the 
HRPT data in 10-bit precision which then are forwarded to the 
NOAA Suitland, Maryland, facility where processing and archiv- 
ing on a more permanent basis occur. The precision loss in 
going from the 10-bit HRPT data to the 8-bit HRPT data is 
between 0.4 and 0.5 degrees when making thermal comparisons 
(Kidwell, 1979) . The amount of data recorded per satellite 
pass depends on the satellite's elevation in relation to the 
receiving station antenna and can be limited by the 13-minute 
capacity of a standard length magnetic tape. Passing directly 
over Redwood City's antenna, NOAA-6 would be within reception 
range for 15.5 minutes, depending on orbital altitude, and 
could provide data from a circular area 6200 kilometers in 
diameter centered on the antenna (Schwalb, 1978) . According 
to Schwalb, the satellite provides useful data only if it is 
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at least five degrees above the horizon. This reduces the 
contact time to 13 minutes and the circular area to 5200 
kilometers . 

a. Field-Station Format 

The field-station format differs from the CDA- 
station format in that it is a combination ASCII-Binary format 
consisting of a single header record at the beginning of the 
tape followed by up to 15000 data records. Each of the data 
records is a sequential interleavening of the scan lines and 
the recorded channels as shown in Table 3 below: 

Table 3 

Field-Station Format 

RECORD CONTENTS 

1 header 

2 scan line 1 — AVHRR 

3 scan line 1 — AVHRR 

4 scan line 1 — AVHRR 

5 scan line 2 — AVHRR 



channel A 
channel B 
channel C 
channel A 



14998 scan line 5000 — AVHRR channel A 

14999 scan line 5000 — AVHRR channel B 

15000 scan line 5000--AVHRR channel C 

channel A, B, or C = any sequence of channels 1, 2, 3, 4 
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The 40-byte header record, all in ASCII, contains 
the ground station identification (SFO for Redwood City) , the 
channel numbers identifying which three of the four available 
AVHRR channels were recorded, the time (GMT) of the first scan 
line, the duration of the pass, and the orbit number. See 
Figure 16 for an example of the header record. Each of the 
remaining 15000 or so data records have identical 2138-byte 
formats beginning with a 14-byte ASCII "mini-header" consisting 
of an identification sequence, the specific AVHRR channel num- 
ber from which the data in the record originated, the Julian 
date of the scan line, and the time (GMT) of the scan line. 
Following the "mini-header" are 10 bytes of telemetry data, 

6 bytes of bac)c scan data, 10 bytes of space view data, and 
50 bytes of space data, all in binary format. The remaining 
2048 bytes, also in binary, are the video data from which esti- 
mates of the sea-surface temperature are derived. See Figure 

17 for an example of one of these data records . 

b. Ephemeris Data Set 

A set of ephemeris data for NOAA-6 also is main- 
tained at Redwood City. An ephemeris data set consists of 
tracking information so that the field station can capture the 
satellite's data stream as the satellite rises above the hori- 
zon and passes overhead to the opposite horizon. See Figure 

18 and Figure 19 for examples of an ephemeris data set. More 
importantly to this project, the ephemeris also contains the 
subsatellite points for the pass calculated at one minute 
intervals. A subsatellite point is that point on the earth's 
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Figure 16 



Header record — field-station format 







(adapted from 


Kidwell , 


1979) 
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BYTE 


CONTENTS 


BYTE 


CONTENTS 


TYPE 


1 


1 


station ID 


2 


station ID 


ASCII 


2 


3 


station ID 


4 


blank 


ASCII 


3 


5 


blank 
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channel A 


ASCII 


4 


7 


channel B 
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channel C 
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hours 
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hours 
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minutes 
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minutes 
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13 


seconds 
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seconds 


ASCII 


8 


15 


duration -min 


16 


duration -min 
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channel A, B, or C 


= channel 1, 2, 3, or 


4 


As an 


example. 


a pass selected 


for the 


project may have a 


header 


record 


as follows: 









SFO 1340333481300 7244 

indicating a Redwood City tape (SFO) containing AVHRR channels 
1, 3, and 4. Time of the first scan line was 03 hours 33 
minutes and 48 seconds (GMT) while the duration of the pass 
recorded was 13 minutes and 00 seconds. The orbit number was 
7244. 
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Figure 17 



Data record — field-station format 
(adapted from Kidwell, 1979) 
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QLUTION NUMBER AT EPOCH TIME .. = 7361 



an example (from Breaker, 1980) 
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surface directly beneath the spacecraft and represents the 
middle of the scan line for the AVHRR. All the subsatellite 
points for each scan line taken together represent the ground 
track that the satellite followed in its orbit. It should be 
noted that any alignment errors made when the AVHRR module was 
attached to the spacecraft during construction may result in 
the subsatellite point not being the center of the AVHRR scan 
line. A summary of alignment data may be found in an undated 
report prepared by ITT Aerospace for NASA. For the purposes 
of this project it was decided that any alignment errors were 
so slight as to be negligible and therefore that the sub- 
satellite point would represent the center of each AVHRR scan 
line. Other information in the ephemeris data set important 
to this project were the orbital elements listed in the preface 
to each ephemeris including orbital period, semi-major axis, 
eccentricity, and inclination. This information was vital to 
the orbital calculations made further on in this project and 
will be explained there. 

c. Initial Satellite Pass Selection 

Each AVHRR scan line is approximately 2840 kilom- 
eters long with 1420 kilometers on each side of the subsatellite 
point. With this information as well as the ephemeris data 
sets for the dates of the P3-C flights and a chart of the 
Northeast Pacific Ocean, it was relatively easy to determine 
specific passes which viewed the ocean areas where the AXBT ' s 
were dropped, thus satisfying the first decision criterion 



64 



mentioned above. Twenty-two satellite passes were thus 
selected for further screening. 

The selection from these 22 passes of orbits whose 
time matched as closely as possible the time of the AXBT drops 
was done in conjunction with the investigation of cloud cover- 
age over the ocean area of interest. As the project relied 
solely on the use of the infrared channels of the AVHRR and 
because cloud cover effectively prevents AVHRR scan coverage 
of the ocean surface, the absence of cloud cover in the ocean 
area of interest was a major factor in pass selection. Also, 
as will be discussed below, NOAA-6 is a sun-synchronous satellite 
that circles the earth 14.2 times every 24 hours. As a result, 
NOAA-6 views the same earth location at the same local sun 
time each day. This translated to our ocean area of interest 
as between 0330 and 0400 (GMT) for ascending passes and between 
1650 and 1720 (GMT) for descending passes. Since Redwood City 
maintains hourly pictures taken from the visual channels of 
the geostationary GOES-WEST satellite, examination of these 
pictures for cloud coverage in the ocean area of interest re- 
sulted in the selection of one pass for each of the six flight 
dates that represented the best compromise between matching 
times and cloud coverage. The six passes chosen for further 
examination are listed in Table 4 below. 

2 . NASA-Ames Research Center 

The last criterion to be satisfied, identification of 
a landmark on each image, was done on the Interactive Digital 
Image Manipulation System (IDIMS) located at the Technology 
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Table 4 



Selected Satellite Passes 



AIRCRAFT 

FLIGHT 


PROJ 

NO. 


AXBT DROP 

TIMES (DTG) 


ORBIT 

NO. 


ORBIT 

TIMES (DTG) 


TYPE 

P.ASS 


15 


Nov 


80 


1 


151905-160049 


7209 


151658-151711 


D 


17 


Nov 


80 


2 


171932-180236 


7244 


180331-180346 


A 


19 


Nov 


80 


3 


191843-192232 


7266 


191712-191716 


D 


01 


Dec 


80 


4 


011826-020031 


7429 


010344-010358 


A 


03 


Dec 


80 


5 


031826-040127 


7465 


031656-031710 


D 


05 


Dec 


80 


6 


051816-060013 


7486 


050355-050408 


A 








where 


D = descending 


and A 


= ascending 





Applications Branch of the Airborne Missions and Applications 
Division under the Director of Astronautics, NASA Ames Re- 
search Center, Moffet Field, California, 
a. IDIMS 

The IDIMS system is a software package that inter- 
acts with a minicomputer (HP-3000), a display terminal, a 25- 
inch COMTAL display screen, a Dunn Instruments color camera 
recorder, and a high-speed printer, and is used extensively 
to work with satellite data, especially LANDSAT imagery. Op- 
tions are available that allow the user to manipulate inter- 
actively satellite imagery so that specific topics of interest 
may be investigated like land use with LANDSAT data, sea surface 
temperature, cloud cover, or ice pack coverage as examples from 
TIROS-NOAA imagery or the many other applications available 
from NIMBUS-7 imagery. The mechanics behind specific options 
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are proprietary data owned by Electromagnetic Systems Labora- 
tory, Inc. (ESL) of Sunnyvale, California, who developed IDIMS 
and to whom the reader is referred for more detailed information 
(ESL, Inc. , 1978) . A second IDIMS system is located at the 
Scripps Institution of Oceanography in La Jolla, California, 
where initial training was done by the author in the use of the 
IDIMS system. A third IDIMS system, also operated by NASA Ames, 
is located in a mobile van that tours the Western United States; 
its facilities were used upon one of its stops at the Naval 
Postgraduate School. 

b. Landmark Identification 

The basic procedure for landmark identification 
used by this project on the IDIMS system was as follows: 

(1) run a short tape routine on each of the se- 
lected pass's magnetic tapes to identify the number of data 
records, hence the number of scan lines per pass (number of 
data records minus 1 header record, all divided by 3) ; 

(2) read the 3-channel AVHRR data from the magnetic 
tape into computer memory and initiate IDIMS processing; 

(3) recall that data comprising the infrared 
channel from memory and display it on the COMTAL display 
screen; 

(4) enhance the displayed infrared imagery for 
temperature using false colors; 

(5) use the Dunn color camera recorder to produce 
an 8 by 10-inch color Polaroid photograph of the enhanced 
image ; 
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(6) use the ZOOM option of IDIMS to enlarge se- 
lected sections of the displayed image in order to locate 
landmark pixels by scan line number and sample number; and 
finally, 

(7) use the PICPRINT option of IDIMS to dump to 
the high-speed printer the count values of all pixels within 
a specified area surrounding the landmark. 

An explanation of certain aspects of this procedure 
is explained in the sections below. 

(1) Count Values. In order to display a satellite 
image on the COMTAL display screen, the IDIMS system sequen- 
tially unpacks the video data read into memory from the 
magnetic tape. These data consist of count values between 0 
and 255 which represent the difference in detected energy be- 
tween a look at deep space and a look at a radiating surface 
such as the earth or clouds (Schwalb, 1978). These count 
values are used by the IDIMS system to produce an image with 
a grey-scale intensity range from 0 to 255 in order to match 
the same range of count values. Options available on the IDIMS 
system allow various color assignments based on these count 
values including an automatic full-spectrum false color assign- 
ment where red is "hot" and blue is "cold" or vice versa. Also 
available are options allowing a single color to vary in its 
saturation over the full count range of 256 values or the 
assignment of a specific color to an individual count as 
"blue=count 125" or to a range of counts as "blue=counts 125 
through 130". Assigning colors in this manner tends to produce 
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a confusing image if many hues are used indiscriminantly . 

For purposes of this project, the automatic full-spectrum 
false color assignment of red "cold" and blue "hot" was used 
so that the oceans were blue and the cloud tops, being much 
colder, were red in all the photographs, 
c. Pixel Identification 

The identification of a landmark pixel and the 
subsequent assignment of a scan line number and sample number 
are made easy by the IDIMS system but the principles behind 
their assignment had to be understood so that other landmarks 
and buoy positions could be located as needed later on in the 
project. 

As seen in Table 4 above, three of the six selected 
NOAA-6 passes were ascending passes and three were descending 
passes. An ascending pass is one where the satellite in its 
orbit crosses the earth's equator heading northwards while a 
descending pass is one where the satellite crosses the equator 
heading southwards. As the satellite is moving, the AVHRR 
scan mirror sweeps from right to left perpendicularly across 
the satellite's velocity vector six times per second with each 
sweep defining a single scan line. Although the scanning 
mirror rotates in a complete circle, only the data located 
55.4 degrees either side of nadir is retained. Nadir is a term 
similar in meaning to the subsatellite point in that nadir 
represents that point on the sweep of the scanning mirror 
when the mirror is pointed at the spot on the earth's surface 
directly underneath the spacecraft. The radiometer data 
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stream from this 110.8-degree arc is electronically divided 
by the MIRP into 2048 equal samples each representing 0.054 
degrees of the total arc (110.8 divided by 2048). Therefore, 
when the satellite transmits its data stream to the earth 
receiving station, the first 2048 bytes of video data taken 
together is intrinsically labeled scan line mamber 1 while the 
first byte is intrinsically labeled sample number 1 or, as 
used throughout this project, (1,1). The second byte of the 
first scan line is labeled (1,2). Sample number 1505 from 
scan line number 3520 would be labeled (3520,1505). Sample 
number 1 through 1024 are located to the right of the sub- 
satellite point when looking in the direction of the satellite's 
velocity vector while samples 1025 through 2048 are to the left. 
Because of this right-to-left pixel numbering system, when an 
image is displayed on the COMTAL unit by the IDIMS system, an 
ascending pass image looks reversed and upside down while a 
descending pass image looks normal where normal is defined as 
having Alaska to the north or top of the image, Hawaii to the 
west or left of the image, aind California to the east or right 
of the image. An ascending pass image, when displayed on a 
display screen or when stored into computer memory for further 
processing, has Alaska on the bottom, Hawaii on the right, and 
California on the left of the image. This occurs because IDIMS 
always displays pixel (1,1) in the upper left corner of the 
COMTAL unit. Although landmark identification can be easily 
done on either type of pass, the geometry involved further on 
in this project rests heavily on a clear understanding of which 
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type of pass you are analyzing and on which side of the sub- 
satellite point is the landmark or buoy. Figure 20 is an 
example of an ascending pass while Figure 21 is an example of 
a descending pass. In each figure, the satellite would travel 
directly up or down the center of the image respectively. 

The IDIMS system automatically keeps track of this 
numbering system and conveniently displays the scan line num- 
ber, sample number, and count value for the pixel you have 
identified on the COMTAL unit using a movable cursor. By use 
of the ZOOM feature on IDIMS and with reference to a chart of 
the local area, it is usually easy to identify landmarks. In 
most of the passes used on this project, the San Francisco Bay 
Area was identified readily and its (scan line number, sample 
number) determined. For purposes described later on, up to 
20 landmarks per satellite pass along the west coast of the 
United States from Glacier Bay, Alaska, south to Mexico and 
from San Francisco east to Pyramid Lake, Nevada, were identi- 
fied. Only one of these landmarks is needed to navigate the 
image as will be explained later on. 

The PICPRINT feature of IDIMS allowed the dumping 
to a high-speed printer of the count values of the pixels sur- 
rounding the landmark. This was done as a method of verifying 
the accurate position of the feature chosen for landmark iden- 
tification. Figure 22 is an example of a PICPRINT output where 
by using a variation of the game of connecting the dots one 
can connect count values in order to recognize features. This 
method of verification will not work, or is made more difficult, 
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Figure 20. NOAA-6 ascending pass, from IDIMS processing 
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Figure 21. NOAA-6 descending pass, from IDIMS processing 
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Figure 22. PICPRINT output of count values, San 

Francisco Bay entrance; cloud-free image 
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if the temperature of the land is the same as the water, if 
any low-lying clouds (fog) have similar temperatures as either 
the water or land, or if clouds obscure sections of the land. 
Normally, the visual channels of the spacecraft are used but 
since many of the passes occurred at night, only the infrared 
channels were usable. Figure 23 is an example of a case 
where clouds interfered with the landmark identification 
process . 

Because an average pass contained 4680 scan lines 
each with 2048 samples or about 9.6 million pixels per pass, 
a convenient method of matching pixel numbers (hence count 
values) to AXBT positions was necessary; thus a system was 
needed to "navigate" the image. 

3 . Satellite Image Navigation 

As mentioned in the introduction, there are many sources 
of error when one wants to compare satellite-derived sea sur- 
face temperatures and AXBT-derived sea surface temperatures. 

It was decided at the beginning of this project that an at- 
tempt would be made to reduce as much as possible one of 
these, that being the earth location errors associated with 
transferring AXBT drop positions to a satellite image, so that 
temperature comparisons could be made. Several methods were 
tried and three of those methods that produced the smallest 
errors are described below. For the purposes of developing 
procedures for locating AXBT's on the satellite image only, 
an assumption was made that the geographical location of the 
AXBT's was accurate; hence any aircraft navigation errors. 
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Figure 23. PICPRINT output of count values, San 
Francisco Bay entrance; cloud-covered 
image. (Note: This image is from an 

ascending pass hence it appears upside 
down . ) 
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ballistic errors on the falling AXBT ' s once launched from the 
aircraft, drift errors on the floating AXBT, or human errors 
in transcribing positional data from aircraft displays to logs 
were ignored. These sources of error will be discussed later, 
a. Zoom Transfer Scope 

A Bausch and Lomb Zoom Transfer Scope was used 
initially in an attempt to transfer the AXBT positions to 
the satellite image. A zoom transfer scope allows one opti- 
cally to overlay a chart, on which the AXBT positions have 
been plotted, onto a satellite image where enough distinguish- 
ing features (landmarks) are evident so that by optically 
stretching, condensing, or rotating the chart, landmarks on 
both chart and image coincide. Once the landmarks coincide, 
the operator manually marks with a pencil the AXBT positions 
onto the satellite image. While the system works fine with 
small area images consisting mostly of land, it could not be 
used satisfactorily on this project for a number of reasons. 
First, each of the selected NOAA-6 passes covered an area ex- 
tending from Northern Mexico to Alaska and from mid-Pacific 
to the western United States. Reducing a chart of this area 
to a size suitable for use on a zoom transfer scope (about 10 
by 10 inches) necessarily requires reduction in the accuracy 
of plotting geographical coordinates. Second, on each of the 
NOAA-6 passes, approximately 80 to 90 percent of the coverage 
area was open ocean with any visible landmass only on the edge 
of the image. Landmasses on the edges of these images are 
much more distorted than landmasses near the subsatellite 
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points due to a combination of satellite scan geometry, earth 
curvature, and the transfer of these images to a flat medium 
like a photograph or chart. Third, the AXBT's were dropped 
along a line over 1300 kilometers long stretching northwestard 
from Cape Mendocino, California. There are no landmarks in 
the Northeast Pacific Ocean between Cape Mendocino and the 
Aleutians; therefore location accuracy decreased the farther 
away from the coast the AXBT's were dropped. Fourth, marking 
a chart manually with a pencil necessarily involves inaccura- 
cies especially when one is trying to locate geographically 
an item as small as an AXBT. Last, and the hardest to over- 
come, is that once the buoy position is marked on the satellite 
image, some method must be found to determine the pixel number 
and hence the count values of the AXBT's position. Remember- 
ing that there are 9.6 million pixels per image, determining 
the exact pixel to choose for a count value would involve some 
guesswork and possibly even large-scale pixel averaging. 
Satellite images unfortunately do not come marked with latitude 
and longitudes, nor do charts contain scan line numbers and 
sample numbers. 

b. IDIM's TRNSFORM 

A second method of trying to locate an AXBT on 
the satellite image involved the use of an ESL, Inc-developed 
IDIMS function called TRNSFORM. TRNSFORM is used mainly in 
registering LANDSAT imagery and involves the calculation of a 
transformation matrix between matching sets of control points 
using a least-squares fit method (ESL, Inc., 1978). A first. 
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second, or third order transformation is possible. TRNSFORM 
was not designed to navigate NOAA-6 imagery mainly because 
TRNSFORM requires 15 to 20 landmarks spread over the entire 
image in order to obtain a small pixel error. Therefore, use 
of this function also failed to provide the accuracy desired 
for this project for one of the same reasons that the zoom 
transfer scope failed, in that landmasses were present only 
on the edges of the images . 

The method finally used, and from which a location 
accuracy of less than 2 pixels resulted, was the development 
of a computer program that determined a satellite's orbit 
referenced to a single landmark in the image and from this, 
when given an AXBT latitude and longitude, could determine the 
scan line nimber and sample number of the AXBT. The development 
of the program required a basic understanding of the orbital 
dynamics of NOAA-6 as well as a working knowledge of spheri- 
cal geometry. 

c. NOAA-6 Orbital Dynamics 

For orbital information in this section, the work 
by Stewart (1979) and Schwalb (1978) was used extensively. 

The NOAA-6 satellite is a sun -synchronous satellite 
which means that its orbital plane rotates at the same rate 
as the rotation of the earth about the sun. As a result, the 
satellite views a point on the Earth's surface at the same 
local sun time each day. Table 4 above listed those times 
that NOAA-6 viewed the AXBT drop area. According to Schwalb 
(1978) , the orbital plane precession rate is approximately 
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equal to 0.000000199 radians per second or 0.986 degrees per 
day eastwards. This rate is achieved by placing the satellite 
in an orbit with a suitable inclination. In the case of the 
NOAA-6 satellite, the inclination was determined prior to 
launch to be 98.739 ± 0.15 degrees, where inclination (i) is 
defined as the angle the satellite's orbital plane makes with 
the earth's equatorial plane measured counterclockwise from 
east. A retrograde inclination (i_) is the supplement of 
the inclination. See Figure 24. 

Figure 24 

NOAA-6 orbital plane inclination 




i = inclination (from ephemeris) 
i = retrograde inclination 

The period of the satellite, obtainable from the 
ephemeris data set (as is the inclination) , is the amount of 
time it takes the satellite to make one orbit of the Earth. 
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The predetermined launch value for the NOAA-6 period was 101.58 
minutes; therefore, NOAA-6 orbits the Earth 14.18 times per 24 
hours. For each orbit the earth rotates 25.40 degrees east- 
ward. See Table 5 for a summary of NOAA-6 orbital parameters. 



Table 5 

NOAA-6 Orbital Parameters 



orbital plane precession rate 
inclination (i) 
retrograde inclination (i_) 
period 

orbits per day 

earth rotation per orbit 

orbital altitude 



0.986 deg/day east 
98.739 ± 0.15 deg 
81.261 ± 0.15 deg 
101.58 minutes 
14.18 

25.40 degrees east 
833 r 18.5 km 



The predetermined launch altitude of the satellite 
was 833 = 18.5 kilometers. The orbit of NOAA-6, to a first 
approximation, is an ellipse. From the ephemeris data set 
the semi -major axis of the ellipse and its eccentricity can 
be found, thus making the satellite's altitude on any pass 
simple to calculate as shown in Figure 25 below. The mean 
satellite altitude (H) is that distance measured from the 
center of the Earth and can be derived using Equation (1) . 

Using the above orbital information, the data in 
Table 6 below was derived for the six selected NOAA-6 passes 
used in this project. 
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Figure 25 

Satellite Altitude Determination 




a = semi-major axis (from ephemeris) in nm 

e = eccentricity (from ephemeris) 

2 1/2 

b = semi-minor axis = a[(l-e )] 

H = mean satellite altitude = (2b+a)/3 (Eq. 1) 



Table 6 

Satellite Data Set Orbital Parameters 



ORBIT 

PARAMETER 



7209 7244 7266 7429 7465 7486 



inclination ( degrees) 


98.69708 


98.69708 


98.69708 


98.69123 


98.69123 


98.69123 


retrograde inclination 


81.30292 


81.30292 


81.30292 


81.30877 


81.30877 


81.30877 


period (minutes) 


101.13285 


101.13285 


101.13285 


101.13084 


101.13084 


101.13084 


semi-major axis (km) 


7185.4875 


7185.4875 


7185.4875 


7199.1856 


7199.1856 


7199.1856 


eccentricity 


0.001187 


0.001187 


0.001187 


0.000603 


0.000603 


0.000603 


mean satellite altitude 


7185.4841 


7185.4841 


7185.4841 


7199.1847 


7199.1847 


7199.1847 
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The reason that the orbital parameters are not 
constant for each pass is that the satellite is subject to 
many forces that tend to cause its orbit to vary. The largest 
of these forces is the fact that the earth is not a perfect 
sphere but an oblate spheroid. King-Hele (1958) and Brouwer 
(1959) developed mathematical solutions to describe this per- 
turbation whose primary effects on the orbit include changing 
the orbital plane precession rate and changing the period. A 
secondary influencing factor is the effect of atmospheric 
drag on the satellite which acts to change the eccentricity 
and is a function of the satellite's altitude. A third influ- 
ence is the effect of solar wind and radiation. It should 
be noted that during 1980, the International Solar Maximum 
Year, solar flare and sunspot activity reached some of the 
highest levels recorded (Ponte, 1981) . Lesser influences in- 
clude the gravitational effects of the sun and the moon on 
the satellite. 

d. Computer Navigation Program 

The main computer navigation program was developed 
under the following premise: given the orbital parameters of 

NOAA-6, the latitude and longitude of an AXBT, and a satellite 
image upon which one landmark has been identified as to (scan 
line number, sample number) and latitude and longitude, calcu- 
late the (scan line number, sample number) of the AXBT so that 
the count value, hence the sea surface temperature, of that 
pixel can be identified readily either on IDIMS or any other 
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computer system. Procedural methods were outlined by Mueller 
(1981) . 

(1) Preliminary Programs. Three preliminary com- 
puter programs were designed to be run on the IBM 3033. The 
first program, SCANLINE, was a simple block counter that counted 
the number of data records on each magnetic tape. The number 

of data records minus the header record divided by 3 gives the 
total number of scan lines per pass. The program listing can 
be found in Appendix B. The second program, TAPEDUMP, was a 
routine designed to dump from the magnetic tape any number of 
bytes per data record and to translate their ASCII-Binary for- 
mats into decimal notation. This program was used to verify 
that there were indeed six scan lines per second and its 
listing can be found in Appendix C. The third program, AREAMAP , 
was designed to function in a manner similar to the IDIMS ' 
PICPRINT function in that it would extract from the magnetic 
tape the count values of a selected grouping of pixels around 
the landmark or AXBT pixel. This program was used to verify 
landmark locations and to determine the surface thermal struc- 
ture around the position of the AXBT. Its listing can be 
found in Appendix D. 

(2) Common Case Geometry. In the development of 
the main computer program, it was necessary to consider four 
cases in the process of predicting an AXBT pixel. These four 
cases are: 

(a) an ascending pass where the landmark 
has a sample number greater than 1024 (Case 1) ; 
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(b) an ascending pass where the landmark 



has a sample number less than or equal to 1024 (Case 2) ; 

(c) a descending pass where the landmark 
has a sample number greater than 1024 (Case 3) ; and 

(d) a descending pass where the landmark 
has a sample number less than or equal to 1024 (Case 4) . 

Common to each of these four cases was the 
assumption of a spherical earth with a radius equal to the 
earth's radius at the landmark latitude. This local radius 
can be calculated using Equation (2) 

2 2 
cos (L ) sin (L ) , 

p — r O j. ^ 1 l/'^ (Pq 

^ ~ M3443.925) (3432.381)^ ^ 



where : 

R = local earth radius in nm 

L = landmark latitude in degrees, 
o 

Also common to all four cases were the calcu- 
lations to determine the great circle distance between the 
subsatellite point and the pixel containing the landmark. 

These calculations refer to Figure 26 below. 

These calculations are made possible under 
the assumptions that the subsatellite point is directly be- 
neath the satellite on the earth's surface, that the scanning 
mirror of the AVHRR forms a scan line perpendicular to the 
satellite's velocity vector, and that the earth is a perfect 
sphere . 
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Figure 26 

Determination of great circle distances 




where; H = mean satellite altitude from Eq. 1 
R = local earth radius from Eq. 2 
9^ = scan angle 
cb = great circle distance 

g 

9g = geocentric angle 

9 = zenith angle 

P 



Determination of the scan angle (9g) degrees 
assumes an equal division of the arc viewed by the radiometer 
(110.8 degrees) into 2048 samples, thus 

((landmark sample-1024) ( 110 . 8/2048) if sample > 1024 \ 

1 (102 5-landmark sample) (110.8/2048) if sample 1024 ) 



The zenith angle (9^) in degrees now can be found: 



(R+H) (sin 9 ) 
1 R 
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The geocentric angle (6^) in degrees is found from Equation 
(3) . 



0 

g 



0 

p 



0 



s 



(Eq. 3) 



If desired, the geocentric angle in degrees can be expressed 
as the great circle distance in nautical miles from Equation 
(4) . 



9 = 60 0 (Eq. 4) 

g g 

(3) Units and Notation. Because the IBM 3033 
and the Fortran computer language were used heavily during 
this project, all angles were converted to or used in radians. 
Table 7 lists the common conversion factors used. Notations 
on all figures included in this project were designed to have 
the same definition whenever possible so comparisons could be 
made between the four cases. 



Table 7 

Program Conversions 



45 degrees = ^ radians = (8) [(tan ^(1.0 radians)] 

, . j. (same angle in degrees) , 

any angle in radians = (tt/4) 



any geocentric angle in 
radians expressed as a _ 45 

great circle distance t/4 

in nautical miles 



(geocentric angle) (60) 



1 nautical mile = 1. 835 Icilometers 
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(4) Nodal and Subsatellite Point Calculations. 



The main program, LOCATE, is divided into two sections. The 
first section calculates the orbital characteristics referenced 
to the previously-identified landmark. The desired output is 
the time and longitude of the ascending or descending node. 
Their calculation is dependent on the four cases enumerated 
above and described in detail below. Once the time and longi- 
tude are known, the second part of the program can proceed to 
calculate the pixel number for an AXBT. To begin, the calcu- 
lation of the ascending or descending node and time follows 
for each of the four cases. 

Case 1. Ascending pass with landmark sample n\amber greater 
than 1024. 

The derivation of orbital characteristics 
referenced to a single landmark in Case 1 made use of Figure 
27 be low . 

At time equal 0, the satellite is directly 
over the subsatellite point that has an unknown latitude 
and longitude (X^) . The only known quantities are that the 
scan line that includes the subsatellite point also includes 
the landmark with known latitude ? longitude (^q) ? 

from IDIMS , a known scan line number (NL) ; and sample number 
(NS) . From the ephemeris data set for this pass, the inclina- 
tion (and hence retrograde inclination (i_) ) a-nd the period 
are known. From Equation (3) or (4) , the great circle dis- 
tance * is known. The goal of this orbital set of calcula- 

g 

tions is to find the latitude and longitude of the subsatellite 
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Figure 27 

Case 1 — orbital characteristics 




point, the longitude of the ascending node the 

time of the ascending node. 

By using similar triangles and the Law of 
Sines, the angle (e) can be determined as follows: 



from triangle I 



sin £ 
sin L 



sin p 



from triangle II 



sin (i_ - e) 
sin p 



sin 



hence , 



= tan 



-1 



sin 1 - 



sin ® 

_ 2 . 

Lsin L 



+ cos 1 
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From triangle I 



, sin L 

-If On 

) = sin [— : ] 

o s in £ 



From triangle II, the distance over which the satellite 
travelled between the ascending node and the subsatellite 
point is 



, cos (p 

-1 r On 

). = cos [ : — J 

t cos (p 



Finally the latitude of the subsatellite point (Lg) can be 
determined using Equation (5) and triangle III 



L 

s 



-1 



sin [ sin (i_) sin ((^^) ] 



(Eq. 5) 



For the moment, the rotation of the Earth is ignored, so the 
change in longitude at the equator between the landmark longi- 
tude (X^) and the fixed ascending node longitude — where 

the "f" indicates a fixed or non-rotating earth derived term) 
can be found by solving triangle I as follows: 



AA 



an 



, cos i)) 

-1 r O, 

cos [ = — J 

cos L 



The change in longitude at the equator between the subsatellite 
point and the fixed ascending node longitude can also be found 
from triangle III 



AX = 



COS 



-COS 



cos L 
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Now, the longitude of the subsatellite point (^g) can be 
calculated easily using Equation (6) 

A = A - (AA^ - AA) . (Eq. 6) 

s o an 



The fixed ascending node longitude is now found by 



f 

an 



A 



o 



AA 



f 

an 



To account for the earth's rotation during the time the 
satellite traveled from the ascending node to the subsatellite 
point, the longitude change due to rotation (AA^) is calculated 
and subtracted from the fixed ascending node longitude as 
follows ; 



At (seconds) = — (period in seconds) , (Eq. 7) 

where At is the time the satellite took to travel between the 
ascending node and the subsatellite point. Continuing, 



A A, 



2 TT 

24 



(1. 002738) (At) 



(Eq. 8: 



where 1.002738 is the sidereal day correction factor. Finally, 



an 



'an - "'r 



(Eq. 9) 
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The time of the ascending node can be found by subtracting the 

At in seconds from the landmark time in seconds. Landmark 

time can be read from the magnetic tape using the tape dump 

program mentioned above to dump the data record containing 

the scan line of the landmark. 

Case 2. Ascending pass with landmark sample number less 

than or equal to 1024. 

The orbital calculations for this case are 

very similar to those of Case 1 and have the same goals. Figure 

28 shows the geometry applicable in this case. As in Case 1, 

the known values are the great circle distance <> from Equa- 

tions 3 or 4 and the landmark's latitude (L ), longitude (A ), 

o o 

and the scan line number and sample number (NL,NS) . 



Figure 28 

Case 2 — orbital characteristics 
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Using spherical triangles I and II, the angle 



(e) can be determined as follows: 



from triangle I 



sin (i- + e) 
sin L 

o 



1 

sin di ' 
' O 



from triangle II 



sin s 
sin 



1 

sin ( 1 )^ ' 



therefore, equating the two equations and solving for e yields 



£ = tan 



-1 



sin 1 



sin L 



sin (|) 



- cos 1 



Continuing with triangle II 






o 



. -1 
sin 



sin 






sin 







and the distance the satellite travels from the ascending 
node longitude to the subsatellite point of the landmark's 
scan line ( 4 )^) also can be found by 



-1 r 

= cos [ 



tan (b 



(sin e) (tan 6 ) 
o 



From spherical triangle III, the latitude of the subsatellite 
point can be found using Equation 10 



L 

s 



. -1 
sin 



(sin i_) (sin (b^) ] 



(Eq. 10) 
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The change in longitude between the longitude of the ascend- 
ing node and the longitude of the subsatellite point now can 
be found from triangle III, ignoring the earth's rotation. 



an 



-1 f 

cos [ 



tan L 



(sin i 



) (tan '1)^) 



The angle (a) from triangle III is 



, sin AA,„ 

. -1, an, 

a = sin [ : : ] 

sin 



which can be used in triangle IV to find the change in longi- 
tude between the subsatellite point and the landmark (AA): 



AA = 



Sin 



(cos ct) (sin ) 

-Ir _g_J 



[■ 



cos L 



The longitude of the subsatellite point (A^) now is found 
using Equation 11 



A = A + AA 

s o 



(Eq. 11) 



The fixed earth ascending node longitude is 



an 



an 



(Eq. 12) 



Using Equations (7) and (8) from Case 1 to determine the 
degrees of longitude through which the earth turns during the 
time the satellite travels from the ascending node to the 
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subsatellite point (Z\,\ ) , the ascending node longitude (X ) 
can be found using Equation (9). The time of the ascending 
node is found exactly as in Case 1. 

Case 3. Descending pass with the landmark sample number 
greater than 1024. 

The case for the descending pass is slightly 
different from the ascending pass cases in that the satellite 
is travelling from north to south. The goals and the known 
factors are identical with the ascending pass cases. Figure 
29 below pertains to the spherical geometry applicable in this 
case. Note that the ascending node becomes the descending 
node (dn) in descending pass calculations. 



Figure 29 



Case 3 — orbital characteristics 



north pole 



CatAliitO 




S 



■\ equator. 

Aq 



ground track 
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The calculations for this case are exactly 
the same as those for Case 1 with some exceptions as noted 
below. The fixed earth change in longitude between the land- 
mark longitude and the descending node longitude can 

be determined from 



AA 



f 

dn 



-1 

cos 



cos 

^cos 




(Eq. 13) 



The change in longitude between the subsatellite point longi- 
tude and the descending pass longitude (AA ) can be found from 
triangle III 



AA = cos 



, cos (f) 

\ 

^COS L 



The longitude of the subsatellite point (A^) can be determined 
from Equation 14 

Ag = Aq + (AA^j^ - AA) , (Eq. 14) 

and the fixed earth descending node longitude now also can be 
determined 





Using Equations (7) and (8) , the degrees of longitude through 
which the earth turns while the satellite travels between the 
subsatellite point and the descending node can be determined 
(aAj.) , and from this the rotating earth descending node 
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longitude can be found using Equation (15) 



The time of the descending node now can be determined by 
adding the time calculated in Equation (7) to the landmark 
time . 

Case 4. Descending pass with landmark sample number less 
than or equal to 1024. 

In this final case, the goals and the known 
quantities are the same as in the other cases described above. 
Figure 30 is used to describe the geometry associated with 
this case. 



Figure 30 

Case 4 — orbital characteristics 

north pole satellite 
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The calculations for this case are exactly 
the same as those for Case 2 with some exceptions as noted 
below. The distance the satellite travels from the subsatellite 
point to the descending node, becomes 



-1 

't = 



The latitude of the subsatellite point, L , now can be found 

s 

using Equation (16) as follows: 



L 

s 



sin ^[{sin i_) (sin <|>. ) ] 



(Eq. 16) 



The fixed-earth change in longitude between the landmark 
longitude and the descending node longitude, riow can 

be found by 



AX 



f 

dn 



cos 



-1 



cos 



cos 




/ 



and the change in longitude between the subsatellite point 
longitude and the descending pass longitude, AX, can be found 
from triangle III 



_ •) cos ']) , 



AX = 



cos 



■cos L 



The longitude for the subsatellite point, X^, is now found 
using Equation 17 



X 



s 



(AX 




) , 



(Eq. 17) 
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and the fixed-earth descending node longitude becomes 



A 



f 

dn 



A 



o 



+ AA 



f 

dn • 



Again, using Equations (7) and (8), earth rotation is con- 
sidered now, with the rotating earth descending node longitude 
described by Equation (18) 



^dn 



A , + AA 

dn r 



(Eq. 18) 



The time of the descending pass is determined by adding the 
time calculated in Equation 7 to the landmark time. 

(5) Buoy Pixel Identification. With the calcu- 
lation of the time and longitude of the ascending or descend- 
ing node, the second part of the program can proceed to calculate 
the (NL,NS) of the AXBT . The previous set of calculations 
referenced the orbital characteristics to the landmark pixel 
(remember pixel = (NL,NS)). The geographical relationship 
between the latitudes and longitudes of the landmark and the 
AXBT are known so the purpose of the remaining part of the 
program is to transform this geographical relationship into 
satellite image coordinates of (NL,NS) . Besides the quanti- 
ties calculated above, the only other known quantity is that 
the AXBT has a unique (NL,NS) . 

As a first guess, any arbitrary scan line can 
be chosen to be the "true" scan line containing the AXBT. One 
of the 2048 samples along the "true" scan line could be the 
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sample number of the AXBT. The aim of the calculations below 
is to prove or disprove, geometrically, that the arbitrary 
scan line is the true AXBT scan line and, once the true scan 
line is selected correctly, to calculate the correct sample 
number. 

The procedure to determine the time that the 
satellite recorded the scan line containing the landmark was 
described above. Since the time of the ascending or descend- 
ing node was calculated in the earlier part of the main pro- 
gram, subtracting the two times describes the satellite flight 
time between the particular node and the subsatellite point of 
the landmark scan line. The assumption was made earlier that 
the satellite's orbit can be considered circular with a mean 
altitude (H) and that the period of the satellite was the 
amount of time it takes the satellite to complete one orbit; 
then the flight time between the node and the landmark scan 
line can be described as a great circle distance in radians by 

= 2- (flight time) , ^ 3 , 

^s period ^ 

Another assumption was made earlier, which was proved by using 
one of the preliminary computer programs described above, that 
the satellite records six scan lines per second. Combining 
these factors, it is simple to determine the difference in 
scan lines between the landmark scan line and the first-guess, 
arbitrarily-chosen scan line; divide by 6 to get the time 
difference between the two scan lines; either add to (ascending 
passes) or subtract (descending passes) this time difference 
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from the flight time; and use Equation (19) to determine the 
distance between the particular node and the subsatellite point 
of the arbitrarily-chosen scan line. The latitude and longi- 
tude of this subsatellite point now can be determined using 
spherical triangles. 

Assuming that the arbitrary scan line is very 
close to the true AXBT scan line, the next step is to find 
the correct sample number along this line. Beginning with 
sample number 1, and then taking every 89th sample (1,90,179, 
...,1959,2048) the distance between the subsatellite point 
and the center of the pixel can be found using Equations (1) 
through (4) for each of the 24 sample numbers. With this 
distance it is possible to calculate the latitude and longitude 
of the center of each of these 24 samples by using one of the 
four cases described below. 

Case 1. Ascending pass with sample number greater than 1024. 

This case uses the geometry of Figure 31. The 
angle (e) can be found easily by 



e 



-1 

cos 



cos 



cos 




(Eq. 20) 



from which 



i = cos^[(cosL)(sine)] . (Eq . 21) 

s s 

The distance, in radians, of can be found by using 
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Figure 31 



1 — pixel determination 




sin L 

'p = sin ^ [— . ; 

’^s sin 1 

s 



(Eq. 22) 



hence the latitude of the possible AXBT position 
be calculated 



L 

P 



sin ^[(sin(i^ -p ) ) (sin i„)] 
s g s 



The angles (y) and 



(AXp) 



now can be determined by 



Y 



sin 



cos 



cos 




/ 
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and 



(sin y) (sin ^ ) 

“p = ' COS - L-^ -l '• 



thus the longitude of the possible AXBT position (A^) is 



determined simply by 



A — A + A A . 

p S p 



Case 2. Ascending pass with sample number less than or 
equal to 1024 . 

Case 2 geometry uses Figure 32. 



Figure 32 

Case 2 — pixel determination 
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The angles (z) and (ig) and the distance (4>g) 
are found using Equations (20) , (21) , and (22) respectively. 

The latitude of the possible AXBT position becomes 



L = sin~^ [sin ( sin(i )] 
p s g s 



The angle (AA^) now is found by 



(sin e) (sin (j> ) 

AX = sin ^ [ — = 2_] 

p cos L 



and therefore the longitude of the possible position is 



Case 3. Descending pass with sample number greater than 
1024. 

Case 3, although using Figure 33 for its 
geometry, uses exactly the same equations as in Case 1 with, 
the exception that the longitude of the possible AXBT position 
is found by 



= X 



- A 



Ai A 



D 
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Figure 33 

Case 3 — pixel determination 




Case 4. Descending pass with sample number less than or 
equal to 1024. 

Case 4, although Figure 34 represents its 
geometry, uses exactly the same equations as in Case 2 with 
the exception that the longitude of the possible AXBT posi- 
tion is found by 



In addition to calculating latitude and 
longitude for each of these 24 samples along the arbitrary 
scan line, it also is necessary to calculate the great circle 
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Figure 34 

Case 4 — pixel determination 




distance (d) between the AXBT geographical position and those 
of the 24 samples using Equations 23 

d = cos ^[(sinL, sinL ) + (cosL. cosL )cos(A (Eg* 23) 

bp op P o 

where : 

L, = true AXBT latitude 
b 

Aj^ = true AXBT longitude . (Eq. 23) 

By taking the sample number that has the 
smallest of these great circle distances, creating a bracket 
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89 samples wide on either side of this center sample, and pro- 
ceeding as before through the appropriate case number for each 
of the 180 samples in this bracket, the sample that has the 
shortest great circle distance between itself and the AXBT 
position can be selected. The reason for the wide bracket 
is to allow for earth rotation and curvature whose effects are 
especially noticeable on the edges of the satellite image. If, 
as before, we assume that we had a scan line very close to 
the true AXBT scan line, a 10 by 10 pixel "square" box is 
created around the sample with the smallest great circle dis- 
tance. In reality, this box is not perfectly square due to 
the curvature of the earth and the motion of the satellite 
during the scan sequence. Those boxes near the subsatellite 
point would be more perfectly "square" than those boxes on the 
edges of the image. 

The scan lines on the top and bottom of the 
box as well as the samples on each of the four corners are 
subjected to the same calculations described above to find 
the latitude and longitude of each of the samples on the four 
corners. The geometry of this box is shown in Figure 35. 

The calculations to find the (NL,NS) of the 
AXBT begins with finding the change in the latitudes and 
longitudes with respect to the changes in the sample and scan 
line numbers. 



5L _ 1,^3 ^1 , ^4 ^ 2 , 

3 (NL) 2 ^ 10 10 ^ 
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Figure 35 
Box geometry 




where L. 

1 



latitude of sample on i-th corner 



A . 
1 



longitude of sample on i-th corner 



L 

o 




latitude and longitude of the sample 
having the smallest great circle dist. 



L ,A = AXBT latitude, longitude 
P P 

L^,L^ = check procedure latitudes 



3L _ 1,^2’^1 

3 (NS) 2 ^ 10 10 ^ 



3A _ l/3’^l , ^4~^2, 

3 (NL) 2 ^ 10 10 ^ 



3A 

3 (NS) 



1 -^ 2’'^1 
2 ^ 10 



+ 



^4 ^3 
10 



] 
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Continuing, two intermediate equations are required 



ii 

1 


3A 


(sample 


A) 


__3A 


(scan 


line 


x) ; 


3 (NS) 


3 (NL) 




3L 


( sample 


A) 


3L 


(scan 


line 


x) . 


3 (NS) 


3 (NL) 



Finally, determination of the (NL,NS) of the AXBT can be 
made 



NL 



A B 

3A/8 (NS) " 3L/3 (NS) 
3A/3(NL) 3L/3 (NL) 

3 A/3 (NS) " 3L/3 (NS) 



NS 



A 

3A/3 (NS) 



r3 A/3 (NL) 
4 a/3 (NS) 



(NL) ] 



Rarely will the arbitrary scan line chosen 
as a first guess be close to the true scan line of the AXBT. 

In this case, after the great circle distances have been calcu- 
lated between the initial 24 samples and the AXBT position, and 
the "square" box has been set up, a check is made to see how 
"small" is the smallest great circle distance. As described 
earlier, the satellite follows a ground track as it travels 
poleward that cuts across meridians of latitude at an angle ser 
up by its inclination- This means that scan lines are not 
oriented east-west along degrees of longitude but are oriented 
northwest-southeast (descending pass) or northeast-southwest 
(ascending pass) crossing many degrees of longitude and lati- 
tude. The check involves comparing the latitude of the AXBT 
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with the latitudes L and L, as shown on Figure 35 above. If 

c d 

the arbitrary scan line is very close to the true AXBT scan 

line, for ascending passes, latitude will be south of the 

AXBT latitude while latitude will be north of the AXBT 

latitude. The reverse is true for descending passes. If 

the arbitrary scan line is far away from the true AXBT scan 

line, both L and L, will be either north or south of the AXBT 
c d 

latitude. In this case, the smallest great circle distance 
is converted to an integral momber of scan lines which are 
added to or subtracted from the first-guess arbitrary scan 
line number depending on whether and were both south or 
north of the AXBT ' s latitude respectively. The jump to a new 
scan line initiates the entire procedure again beginning with 
the steps necessary to calculate Equation (19) . This jump 
process terminates when the number of scan lines to be jumped 
is 5 or less at which time the boxing procedure begins with 
the eventual determination of the AXBT's (NL,NS) . 

The main program, whose listing may be found 
in Appendix E, was initially set up to be run interactively on 
a display terminal. The program was used to locate all the 
AXBT's that were dropped from the P-3C. Verification of the 
accuracy of the program was done by selecting 15 to 20 land- 
marks per image and asking the program to predict the (NL,NS) 
even though they were known already from the IDIMS system. 
Results of this verification will be discussed below in Section 
IV. A. 
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C. GOSSTCOMP 



The GOSSTCOMP sea surface temperature charts were obtained 
for the period of this project from NOAA-NESS. These charts 
are produced on a weekly basis by NOAA-NESS using procedures 
outlined by Brower a^. , (1976) and since updated to take 

advantage of the AVHRR on NOAA-6. 
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IV. RESULTS 



A. NAVIGATION ACCURACY 

A major effort was made on this project in an attempt to 
reduce the effects of geometric distortion associated with 
locating landmarks or open-ocean positions on satellite imagery. 
Previously published works with earth location errors in ex- 
cess of 10 kilometers at nadir were suitable for regional 
location and analysis of mesoscale features; however, it was 
believed that accurate comparisons of thermal data were signi- 
ficant only if the products being compared were co-located in 
the same geographical position. The location of thermal 
features is especially important in naval tactical applications. 
If a submarine were taking advantage of the unique acoustical 
properties of an eddy or ocean front, then acoustical prosecu- 
tion by the opposing forces would be more successful if the 
sensors employed by this group were located so as to take 
advantage of the thermal feature also. If pre-mission informa- 
tion mis located the edge of the front or eddy due to the 
geometric distortion inherent in satellite imagery, then the 
results could be disastrous to one of the parties. Ultimately, 
any error in sensor placement could prove disastrous whether 
caused by satellite imagery or not; part of the purpose of 
this project was to make the location error as small as possible. 

Using the main computer program, the (NL,NS) of each AXBT 
was predicted. The error in this prediction was determined 
to be within 2 pixels of the true AXBT (NL,N3) . The procedure 
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to verify this accuracy began with the identification of the 
(NL,NS) of up to 20 landmarks on each satellite image. Each 
of the landmarks then was treated like an AXBT and its geo- 
graphical coordinates were input to the computer to see what 
the program would predict for each landmark's (NL,NS) , These 
predictions then were compared to the IDIMS-determined 
(NL,NS) , and the separations in pixels were determined. The 
results are s\ammarized in Table 8 below. 



Table 8 

Statistical Summary of Navigation Accuracy 



SCAN LINE ERROR 
mean 

standard deviation 
99% confidence level 

SAMPLE ERROR 
mean 

standard deviation 
99% confidence level 



1.32 scan lines 
0.93 scan lines 
0.88-1.90 scan lines 

1.35 samples 
1.17 samples 
0.81-1.90 samples 



The conclusion drawn from this statistical summary was 
that the predicted AXBT (NL,NS) is within 2 pixels of the true 
AXBT (NL,NS). Because of earth curvature, pixels close to 
nadir are not as wide as those out on the edges of the image. 
Sample number 1 and 2048, on the right and left edge of the 
image respectively, are 4.3 kilometers wide whereas sample 
number 1024 and 1025, located to the right and left of nadir 
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respectively, are only 0.77 kilometers wide. As a result, 
the 2-pixel error can be as small as 1.9 kilometers, if the 
predicted (NL,NS) is at nadir; or as large as 10.7 kilometers, 
if the predicted (NL,NS) is at the edge of the image. Table 
9 lists the navigation error associated with selected sample 
numbers. Notice that the error is not linear with distance 
from nadir but is less than 5 kilometers over 80% of the 
image and less than 3 kilometers over 50% of the image. Only 
on the outer 10% of the image does the error balloon from 4.8 
to 10.7 kilometers. 



Table 9 

Navigation Errors Associated with a 2-Pixel Error 



AMPLE NUMBER 


ERROR (km) 


1 


10.7 


200 


4.8 


400 


3.0 


600 


2.5 


800 


2.0 


1024 


1.9 


1200 


2.0 


1400 


2.3 


1600 


3.0 


1800 


4.3 


2048 


10.7 
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There are some other sources of navigation error that 
should be kept in mind when using data developed by this method. 
Alignment of the AVHRR module during its assembly prior to 
launch could be the source of constant offset error. This 
type of error was described previously in this paper; analy- 
sis of the navigation results did not show any consistent off- 
set bias that could be attributed to module alignment errors . 

Through the process of verifying the 2-pixel accuracy, 
many landmarks were identified on the IDIMS system as explained 
above. The data from NOAA-6 infrared channel number 4 were 
used for landmark identification and their use could introduce 
errors in assignment of the (NL,NS) . These errors arise from 
trying to identify ' landmarks whose surface temperature may 
not be very different from the surrounding surfaces. This 
effect would become even more pronounced if ground fog were 
present. When selecting these landmarks, the best contrast 
was effected by land-water boundaries, examples of which were 
Point Lobos west of Monterey, the San Francisco Bay entrance, 
Alcatraz Island, Point Reyes, the Columbia River mouth. Lake 
Tahoe, and Glacier Bay among others. Many possible landmarks 
were not considered if there were insufficient contrast to 
identify the feature. A good example of this was the Seattle- 
Tacoma -Olympia area where the numerous bays and tributaries 
had surface temperatures close to land temperatures , thus 
making it very difficult to distinguish a specific pixel as 
being some peninsula or promontory. 
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A third source of error involved the number of signifi- 
cant figvires used in the mathematics of this project. Single- 
bit precision was used during computer processing. Although 
this may have had an effect after numerous computations (the 
average program run to calculate 24 AXBT positions per image 
executed 475,000 statements), it was felt that the number of 
significant figures was more critical. An example of this was 
the determination of decimal geographical coordinates. The 
navigation system on the P-3C supplied the computer-calculated 
coordinates of the AXBT's to seconds of latitude or longitude. 
One second of latitude error is equal to 0.1 kilometers, 
which in itself is not so large; however, most landmarks were 
identified using charts with scales of 1:2,000,000. After 
determination of the coordinates, a decimal conversion to three 
decimal places was completed. If errors in this procedure 
were compounded by weak land-water contrast on the infrared 
image used for landmark identification, it could contribute 
significantly to the 2-pixel error. 

A fourth source of error has to do with the resolution of 
the AVHRR itself. As discussed earlier, the 1.1 kilometer 
resolution would necessarily make it difficult to identify 
something like the Transamerica Building in downtown San 
Francisco. An example of where this could be a contributing 
factor to the 2-pixel error would be in using the most western 
point of Point Reyes. If the scan sequence is such that the 
radiometer does not resolve this point, then the first pixel 
it does identify as being land would be to the east of the 
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point. The user of the satellite image would have a very 
difficult time in trying to determine whether or not this has 
occurred. As a result, the user would assign the geographi- 
cal coordinates of the most western point, introducing a 1- 
pixel error immediately before any computer processing begins. 

It is felt that if landmark's could be identified using sharper 
land-water contrast or using a visible channel vice the infra- 
red channel if it is available, and if geographical coordinates 
could be assigned with greater accuracy, the majority of the 
2-pixel bias would be eliminated. 

Another error to be considered is that the satellite may 
not be perfectly stable in its orbit. It is likely that small 
amounts of pitch, roll, or yaw occur from time to time although 
the ADACS system was designed to keep these attitudes to a 
minimum. 

The last error to be discussed is the use of several assump- 
tions made during this project. The Earth was assumed to be 
spherical and although the radius was calculated to be that 
radius at the latitude of the landmark, a small error will be 
introduced in calculations involving the earth radius term at 
the latitude of the buoys. Similar small errors arise with 
the assumptions made concerning the satellite orbit during 
scan line calculations, and with the calculation of the mean 
altitude of the satellite above the earth's surface. 

In conclusion, it is believed that the 2-pixel error found 
on this project could be reduced further to sub-pixel accura- 
cies if some of the errors described above were eliminated 
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or refined. Since the development of LOCATE, a program using 
most of the same techniques as LOCATE has been developed to 
predict the geographical coordinates of open ocean images 
with similar accuracies (Mueller, 1981) . 

B. THERMAL COMPARISONS 

1 . Horizontal Distribution 

As usually can be expected when satellite-derived sea 
surface temperatures that are xmcorrected for atmospheric 
attenuation are compared with AXBT values of sea surface 
temperature, the satellite temperatures were colder than the 
AXBT data by a mean difference of 2.9 degrees C. Table 10 
lists the mean temperature difference values and the corres- 
ponding standard deviations for this and following comparisons. 
Figures 36, 37, 38, and 39 show sea surface temperature com- 
parisons of this and other methods to be described below along 
the buoy line. The majority of the 2.9-degree error can be 
attributed to the effects of the intervening atmosphere be- 
tween the ocean's surface and the satellite radiometer. Cloud 
contamination of the satellite values was not considered a 
major factor due to the screening process that went into se- 
lecting the data. Out of the six satellite passes selected 
for study at the beginning of this project, only three met the 
full requirements that were required for processing. Two of 
the passes were not considered due to scattered clouds over 
enough buoy positions to make any comparisons useless and one 
pass was not considered because it contained no clearly 
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Table 10 



Temperature Comparison Statistics 



COMPARISONS 



satellite vs. AXBT 


MEAN (C) 


STANDARD 

DEVIATION 


17 November 


-3.0 


0.5 


01 December 


-2.6 


0.4 


05 December 


-3.0 


0.6 


overall 


-2.9 


0.5 


satellite vs. GOSSTCOMP 


17 November 


-2.0 


0.9 


01 December 


-3.2 


0.6 


05 December 


-4.3 


0.9 


overall 


-3.2 


1.1 



GOSSTCOMP vs . AXBT 



17 November 


-0.7 


0.7 


01 December 


LO 

O 


0.6 


05 December 


1.2 


0.6 


overall 


0.3 


1.0 
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Figure 36. Sea surface temperature comparisons, 
17 November 1980, center track 
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Figure 37. 



Sea surface temperature comparisons, 
1 December 1980, center track 
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Figure 38. Sea surface temperature comparisons, 
5 December 1980, center track 
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Figure 39. Sea surface temperature comparisons, 
5 December 1980, northern track 
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identifiable landmarks. Passes selected for complete pro- 
cessing were 17 November, 1 and 5 December 1980. From personal 
observations onboard the P-3C during AXBT deployment, it was 
noted also that there was no ground fog to interfere with 
satellite measurements of the AXBT positions. The time 
difference of one, three, and three hours between the last 
buoy drop and the satellite flyover for 17 November, 1 and 
5 December, respectively, probably is not a factor as the 
lowest and highest mean error values were found on 1 and 5 
December with 17 November having an intermediate value. If 
there were a correlation, one would expect 17 November to have 
the smallest error but this was not the case. The transient 
warming of the surface waters during the afternoon, the so- 
called afternoon effect, did not occur due to the weather 
conditions during the three-week project period; hence, this 
process also was ruled out as a source of the error. Con- 
stant wind speeds in excess of 20 knots from the south on 17 
November, in excess of 25 knots from the northeast on 1 Decem- 
ber, and in excess of 15 knots from the westnorthwest on 5 
December (National Weather Service, 1980a) along the buoy 
line kept the surface waters under constant wind-mixing condi- 
tions. In addition to the winds, ship observations of the 
sea state at the northern end of the buoy pattern found four 
to ten foot swell and two to six foot waves. These turbulent 
mixing conditions are diametrically opposed to the formation 
of the afternoon effect (James, 1966). 
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A method of "field-calibrating" the satellite data to 
eliminate the effects of the atmosphere was suggested by 
Tabata and Gower (1980) . Using a simple linear regression 
technique, they plotted numerous ship-obtained surface tempera- 
tures versus satellite count values and found that over a 
limited area and a limited time period between satellite and 
ship observations (1.5 days), the error between satellite and 
ship values could be reduced to 0.5 degrees C. This tech- 
nique was tried using the data from this project. An absolute 
mean difference of 0.3 degrees (s = 0.2 degrees) was found 
between satellite and AXBT values. The time period between 
satellite and AXBT observations was three to eight hours. 

Flight crews usually will not have the luxury of ex- 
pending 24 AXBT's on a tactical mission however, so the linear 
regression technique was tried using only two buoys. The 
rationale behind using two buoys was that this is the number 
of AXBT's usually carried on both S-3A and P-3 aircraft. Addi- 
tionally, a scenario could exist whereby a satellite photo 
obtained prior to the flight could pinpoint two sections of 
the tactical operating area where thermal differences exist 
and those two locations could be designated for AXBT deployment. 
The point is to try and get a spread in temperature between 
the two AXBT's. Using the two-buoy method, the data from AXBT 
positions 1 and 13 on 1 December were used for the linear re- 
gression. Predictions of temperatures from count values found 
a mean difference error of 0.3 degrees (s = 0.3 degrees), the 
same value found using the 2 4 -buoy method. When this same 
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regression formula with constants calculated from 1 December 
data was used to predict 5 December (4 days later) and 17 
November (15 days earlier) temperatures, mean errors of 0.5 
degrees (s = 0.3 degrees) and 0.4 degrees (s = 0.3 degrees), 
respectively, resulted. Although examination of all possible 
cases would be necessary before conclusive results could be 
stated, these preliminary estimates indicate that it should 
be possible to use two AXBT's and an infrared satellite image 
to predict sea surface temperatures within 0.5 degrees for at 
least two days after the original satellite pass, a prediction 
tool particularly helpful if clouds obscure the sea surface 
during those two days or if AXBT assets are in short supply. 

In addition, these procedures can be used in near real-time 
processing of current satellite images and do not rely on any 
atmospheric model processing. In any case, these thermal 
predictions that are very accurate in location and fairly 
accurate in temperature would be tactically significant, in 
place of mean values or best-guess values, when doing sound 
velocity calculations near meander, eddy, or frontal regions. 

Relative temperature gradient analysis displayed the 
expected correlation between satellite and AXBT values. Both 
the AXBT and satellite gradients were 0.6 degrees per 60 nm 
on 17 November and 1 December while on 5 December the AXBT 
gradient was 0.56 degrees per 60 nm and the satellite gradient 
was 0.52 degrees per 60 nm. 

An interesting thermal feature whose horizontal sur- 
face manifestation was detected by both the satellite and 
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the AXBT was a meander in the final stages of closing off 
from its parent body of water to form an eddy. See the darker 
region transected by the buoy line through buoy positions 23, 

3, 22, and 4 in Figure 40. This warm-core meander had an 
approximate 100 nm diameter with the exception of an open arm 
extending southward into its parental water mass. The diameter 
was verified by the buoys dropped on the northern and southern 
tracks. These buoys were 60 nm away from the center rack and 
the warm meander did not show up on any of the thermal traces. 
The satellite indication of this diameter resulted in a slightly 
larger radius, a fact attributed to the thermal resolution limi- 
tations of the satellite data in determining weak temperature 
boundaries. In the satellite images, this meander is surrounded 
on the west, north, and east by the Subarctic Current-Calif ornia 
Current confluence. The center of the meander had a surface 
temperature of 15.8, 14.9, and 14.4 degrees C on 17 November, 

1 and 5 December, respectively. A chart of the monthly mean 
surface temperature for November 1980 (Renner, 1981) clearly 
shows the intrusion of a large tongue of warm water from the 
area between San Francisco and Hawaii northward along the 
west coast of the United States. 

The decrease in the surface temperature of the meander 
over the project time period was reflected by the decrease 
in the surface temperatures all along the buoy line. Both 
the satellite and the AXBT ' s recorded mean changes of 0.9 
degrees between 17 November and 1 December and 0.7 degrees 
between 1 December and 5 December. This drop in temperature 
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is to be expected considering the weather conditions as men- 
tioned above. On 17 November, a low pressure system was 
firmly entrenched over the Aleutians while a high pressure 
system was anchored off of Southern California. This pressure 
pattern is typical of the Northeast Pacific in early winter. 

A cold front extending southward from the Aleutian low moved 
across the buoy line during the evening of 17 November and 
crossed the U.S. coastline during the morning of 18 November. 
Winds before passage were southerly at 25 knots while after 
passage the wind shifted northerly at 30 knots; hence the 
condition for considerable wind-mixing existed. A series of 
cold fronts on 20-22 November and 25 November also passed 
through the project area continuing to lower the sea surface 
temperature and drive the mixed layer depth deeper. On 1 
December, low pressure cells were established west of the 
coast of Washington and about 500 nm west of the central 
California coast. The Washington low strengthened and cen- 
tered near buoy positions 7 on 2 December. This low was 
accompanied by winds in excess of 35 knots on 3 December over 
the entire buoy area while a slow-moving cold front hugged 
the coastline. On 4 December, a high pressure area, previously 
established in the Gulf of Alaska, moved into the project 
area from the north pushing the cold front well inland although 
the low remained off the Washington coast. The high moved 
southerly on 5 December, influencing the weather over the 
entire buoy area (National Weather Service, 1980b) . See 
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Figures 41, 42, and 43 for the surface weather depiction 
charts for 17 November, 1 and 5 December, respectively. 

Because the linear regression model mentioned earlier 
was used to remove the effects of the atmosphere with fairly 
good results, a comparison was made between the satellite 
data and the GOSSTCOMP product. Sea surface temperature pro- 
ducts from GOSSTCOMP have been subjected to an atmospheric 
correction model and are issued on a weekly basis. On all 
three days, the satellite data from this project were colder 
than the GOSSTCOMP values. The mean difference for 17 November, 
1 and 5 December were 2.0 degrees (s = 0.9 degrees), 3.2 
degrees (s = 0.6 degrees), and 4.3 degrees (s = 0.9 degrees) 
respectively with the overall mean of 3.2 degrees (s = 1.1 
degrees). This overall mean agrees fairly well with the 3.5 
to 3.9 degree bias enumerated by Klein (1979). The reason for 
the 3.2 degree bias can be attributed directly to the effects 
of the atmosphere, exactly the same situation as seen in the 
AXBT versus satellite comparisons mentioned previously. An 
interesting point to be made is that the 3.2 degree bias of 
GOSSTCOMP versus satellite data is higher than the 2.9 degree 
bias of AXBT versus satellite data. This led to a comparison 
between GOSSTCOMP and AXBT data with the result that GOSSTCOMP 
values were 0.3 degrees (s = 1.04 degrees) warmer overall than 
the AXBT values. Because the project area was never totally 
cloud-free during the period of observations, it is felt that 
the overcorrection for atmospheric effects described by Klein 
(1979) is still a factor in the warmer GOSSTCOMP values. It 
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(from Reed and Mullen, 1981) 
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Surface Analysis for 17 NOV 80, OOOOZ. 



Figure 42. Surface weather chart, 1 December 1980 
(from Reed and Mullen, 1981) 
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Surface Analysis for 1 DEC 80, 00002. 



Figure 43. Surface weather chart, b December ±y«u 
(from Reed and Mullen, 1981) 
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Surface Analysis for 5 DEC 80, OOOOZ. 



should be mentioned that the GOSSTCOMP product did not 
indicate the warm meander that the AXBT and satellite data 
located, probably due to the large grid structure used by 
GOSSTCOMP. 

All three methods of sea surface temperature determin- 
ation, AXBT, satellite, and GOSSTCOMP, were compared to the 
20-year mean surface temperature values of Robinson (1976) . 

AXBT values versus climatology resulted in the November AXBT ' s 
being 0.3 degrees colder than the mean while the December 
values were about the same as the mean. The reason for the 
cooler surface waters in November is probably a result of 
the high incidence of weather frontal passage with accompany- 
ing high winds through the project area. As was expected, 
climatology did not show the warm meander. 

Comparison of satellite versus monthly mean data 
found the satellite data averaging 3.0 to 2.7 degrees colder 
than the mean for November and December. Comparison of GOSSTCOMP 
versus the mean resulted in GOSSTCOMP being 0.4 degrees warmer 
than the mean for November and 2 . 0 degrees colder than the 
mean for December. 

2 . Vertical Distribution 

There is no known way at present to sense remotely 
the vertical thermal structure in the ocean; however, if one 
combines knowledge of the horizontal gradients with clima- 
tology, a fairly accurate synopsis of the upper ocean thermal 
structure is possible. A more accurate picture can be 



134 






^ 



f 











formulated if the satellite data are augmented with well- 
placed AXBT drops. 

From climatology, the expected mean surface tempera- 
ture and the mean layer depth for November were 12.9 degrees 
and 50 meters (s = 5 meters) respectively over the project 
area. The AXBT mean surface temperature and mean layer depth 
for 17 November were 12.6 degrees and 58 meters (s = 6 meters). 
For December, climatology means were 10.8 degrees and 67 meters 
(s = 6 meters) and the AXBT means were 10.8 degrees and 71 
meters (s = 7 meters). Figures 44, 45, 46, and 47 show the 
vertical structure along the buoy line on 17 November, 1 and 
5 December (center track and north track) respectively. From 
the numerous oceanographic studies in the area (Tully, 1961; 
Tabata, 1961; etc.), it is known that during this period of 
the year, the layer depth is deepening towards the maximum 
limit of the top of the permanent halocline at 100 meters. 

The 100-meter depth is not reached usually until February. 

The deepening of the layer is directly attributable to the 
turbulent mixing conditions caused by the sustained high wind 
speeds and by the convective mixing caused by the surface 
cooling during the calmer periods. The weather pattern for 
late-November and ear ly-December was discussed previously. 

A general rule of thumb is that warm surface waters generally 
exhibit shallow layer depths while colder surface waters ex- 
hibit deeper layer depths. This pattern held true throughout 
the project area. Although definitive sea surface temperature- 
layer depth relationships were not within the scope of this 
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Figure 44. Vertical thermal structure, 17 November 
1980, center track 
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Figure 46. Vertical thermal structure, 5 December 
1980, center track 
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Figure 47. Vertical thermal structure, 5 December 
1980, northern track 
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project, it is conceivable that flight crews, after observing 
a horizontal and vertical thermal analysis done on AXBT data 
and satellite images from a previous day(s) could view cur- 
rent satellite infrared images and through comparisons of 
surface temperatures (atmospherically corrected or not) could 
estimate the layer depth. 

Further deployment of layer depth prediction tech- 
niques may be found in combining the works of James (1966) , 
McAlister and McLeish (1970) , and Mollo-Christensen and 
Mascarenhas (1979) . James described a way of using heat bud- 
get calculations and wind-mixing parameters to predict the 
ocean thermal structure. McAlister and McLeish described an 
airborne system that was capable of measurement of the total 
heat flow from the sea. Mollo-Christiansen and Mascarenhas 
used LANDSAT data to calculate heat storage in the upper mixed 
layer of the ocean. The use of satellites to estimate wind 
speed and direction has been demonstrated. If the heat budget 
could be estimated using the principles described in the papers 
above, and the winds determined, then the procedures described 
by James may be applicable. 

An example of how present satellites are inadequate to 
sense the vertical structure is the fact that only the AXBT 
traces located a region of a subsurface temperature maximum 
on the northern end of the buoy line. This region would help 
to define the lower extent of a subsurface sound channel. 

There was no surface manifestation of this feature, like the 
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warm meander, and hence the satellite completely missed it. 
This temperature maximum was relatively narrow (50 meters) 
and had an axis between 150 and 175 meters. On 17 November, 
the area affected by the temperature maximum was relatively 
large with the axis at the shallower depth. By 5 December, 
the area affected was reduced by half with the axis depth 
deepening in conjunction with the layer depth. The area of 
the subsurface temperature maximum can be seen in Figures 44, 
4 5, 46, and 4 7 . 
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V. CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER STUDY 



It was determined that geographical locations over open 
ocean areas could be located on satellite imagery to within 
2 pixels or within 3 kilometers over 50% of the image (within 
5 kilometers over 80%, worst case within 10.7 kilometers). 
Location of satellite features on geographical charts using 
a variation of the main program has similar accuracies. 

Satellite observations can be a very effective tool if 
used in conjunction with available groundtruth data. The 
development of the NOAA-6 AVHRR allows a more accurate de- 
termination of sea surface temperatures than from previous 
satellites; however, the satellite mean values differ from 
groundtruth values by 2.9 degrees C. This bias is attribu- 
table directly to the effects of the atmosphere. By using a 
method of linear regression, it is possible to "field-calibrate” 
the satellite data so that the mean error between satellite 
and AXBT data is reduced to 0.3 degrees. This small bias 
held true whether 24 or 2 buoys were used to define the cali- 
bration constants for the linear regression equation for a 
single satellite pass. When these same constants were used 
to predict the sea surface temperature on satellite passes 4 
days later and 15 days previous, the mean error increased 
slightly to between 0.4 and 0.5 degrees. This three-week pre- 
diction period saw the passage of numerous cold fronts and 
experienced long periods of turbulent mixing conditions so the 
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relatively small bias may make this procedure a possible 
prediction tool in calibrating satellite images without the 
need to use atmospheric attenuation models. 

Relative temperature gradients were constant, as ex- 
pected, between satellite and AXBT data values. The intro- 
duction into the project area of a warm-core meander in the 
final stages of becoming an eddy was apparent readily from 
both the AXBT data and the satellite imagery, although it was 
ignored by GOSSTCOMP . The detection of subsurface thermal 
structures by satellites is successful only if there is a 
surface manifestation of the feature such as the meander. A 
subsurface temperature maximum that is associated with a sub- 
surface sound channel was not detected by the satellite, as no 
surface manifestations were present. 

When GOSSTCOMP values were compared to AXBT values , an 
overall 0.3 degree bias was noted, although individual daily 
values ranged from 0.7 degrees colder to 1.2 degrees warmer 
than the corresponding AXBT data. 

Over the project period, the sea surface temperature 
decreased on the average 1.8 degrees while the layer depth 
deepened on the average 13 meters. The cooling and deepening 
process was directly related to the number of frontal passages 
and the high wind speeds during the project period. It may 
be possible to use systems currently under development to sense 
remotely the existing thermal structure to 100 meters; or to 
use a combination of satellite instruments to determine the 
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heat budget and wind speeds and from this information to 
calculate the thermal structure. 

Determining the exact location of thermal features on 
satellite imagery is only the beginning in tactically employ- 
ing satellite data for naval missions. Further development 
of the relationship between the sea surface temperature and 
the vertical thermal structure is warranted. An excellent 
starting point would be along lines similar to those developed 
under the ASWEPS program. If a possible relationship could 
be derived, the use of environmental satellites by naval 
tacticians could advance far beyond its present stage. 

Using methods developed by this thesis, it is possible to 
develop hypothetical scenarios that need to be tested opera- 
tionally. Supposing that the 17 November satellite imagery 
and AXBT data, and the 1 December satellite imagery only were 
available, an antisubmarine aircraft flight crew assigned a 
surveillance mission on 1 December could have predicted within 
one half degree the sea surface temperature over a wide ocean 
area. Knowing both the vertical thermal structure from 17 
November and the effects of cooler surfaces and sustained 
high winds on this structure, the flight crew could make a 
fairly accurate prediction of the 1 December oceanographic 
and acoustic conditions. Development of techniques that relate 
the surface temperature to the vertical structure would make 
this process that much easier. 

In addition to predicting the vertical structure, some 
aspects of the horizontal structure are equally important. 
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especially to the fast-paced antisubmarine warfare efforts 
of the carrier-based S-3A aircraft. Submarines may be able 
to use sharply-delineated fronts and eddies to keep themselves 
acoustically hidden from aircraft carriers while remaining 
within weapons firing range. The persistence of these thermal 
features over a period of time is seen easily in satellite 
images and their exact geographical location can be determined 
using the methods derived in this thesis. A current satellite 
photo is far superior to other products now available in con- 
veying this type of information. As an example, GOSSTCOMP 
missed totally the small (100 nm diameter) warm-core meander 
found in this project. GOSSTCOMP is also not a real-time 
product . 

In conclusion, the following recommendations for further 
study are suggested: 

(1) the development, where applicable, of an empirical 
relationship between the sea surface temperature and the 
vertical thermal structure; 

(2) the development of thermal structure prediction tech- 
niques using both satellite data and the empirical relation- 
ship developed above; 

(3) a study of the persistence of horizontal thermal 
features using satellite imagery; 

(4) a study of how accurately satellite-observed surface 
thermal features reflect the subsurface structure; 

(5) the development of a streamlined LOCATE program 
suitable for ship-board use so that surface thermal reatures 
from satellite imagery can be located more accurately; 
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(6) the further development and testing of the "field- 
calibration” technique of dealing with atmospheric attenua- 
tion; and, 

(7) the investigation of including a current satellite 
image, with thermal features geographically located, in 
mission planning information. 
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APPENDIX A 



COUNT- 


-TO -TEMPERATURE CONVERSION TABLE 
(from Kidwell, 1979) 


COUNT 


TEMPERATURE (degrees C) 


95 


16.33 


96 


15.91 


97 


15.48 


98 


15.06 


99 


14.63 


100 


14.20 


101 


13.77 


102 


13.34 


103 


12.90 


104 


12.03 


105 


11.59 


106 


11.59 


107 


11.15 


108 


10.70 


109 


10.26 


110 


9.81 


111 


9.36 


112 


8.91 


113 


8.46 


114 


8.00 


115 


7.54 
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116 

117 

118 

119 

120 

121 

122 

123 

124 

12 5 



7.08 



6.62 

6.16 

5.69 

5.23 

4.76 

4.28 

3.81 

3.33 

2.85 
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WRITE(6f300) TAP01340 
URITE<6r301) TAP01350 
CONTINUE • TAP01360 
STOP TAP01370 
END TAP01380 
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LOGICAL*! A(l)rB(4) ARE01400 
B<4)=A(1) ARE01410 
RETURN ARE01420 
END ARE01430 

ARE01440 



C***** JOB CONTROL STATEMENTS ***** ARE01450 
C***** TAPENAME IS THE NAME OF THE MAGNETIC TAPE ASSIGNED WHEN THE ARE01460 
C***** TAPE IS CHECKED INTO THE COMPUTER CENTER. ARE01470 
C ARE01480 
//GO.METTAP DD UNIT=3400-3 r VOL=SER=TAPENAME r DISP=OLDr LABEL=( 1 » BLP) r ARE01490 
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